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1.  INTRODUCTION 


Through  the  decade  of  the  sixties,  computation  of  radiation  environ- 
ments from  nuclear  weapons  was  limited  by  computer  technology  and  the 
paucity  of  basic  data  in  the  form  of  reaction  cross  sections.  These  limit- 
ations resulted  in  high  degrees  of  uncertainty,  i.e.  factors  of  3 and  more, 
in  the  results  of  all  radiation  transport  calculations,  including  those  for 
weapons  radiation  output  and  subsequent  transport.  In  fact  confidence  in 
calculated  weapons  output  was  so  low  that  experimentally  determined 
values  were  preferred  as  sources  in  weapons  radiation  environment  problems, 
while  in  tne  case  of  radiation  transport  in  the  atmosphere  gross  experiments 
such  as  HENRE^1  were  attempted  in  order  to  gain  confidence  in  the  then 
state-of-the-art  of  transport  technology.  The  lack  of  success  in  instilling 
confidence  in  results  was  due  not  only  to  limitations  in  calculational  cap- 
abilities but  also  to  the  ambitious  nature  of  these  experiments,  which, 
when  combined  with  the  then  state-of-the-art  in  experimental  technique, 
produced  measured  results  which  were  difficult  to  reproduce  or  interpret 
in  the  context  of  radiation  transport  calculations.  The  rather  large  un- 
certainties in  radiation  environments  led  the  Defense  Nuclear  Agency  (DNA) 
to  support  two  major  technology  programs,  one  to  improve  the  calculation 
of  weapon  radiation  output  and  the  other  to  improve  radiation  transport  in 
the  environment,  with  special  emphasis  on  transport  in  the  atmosphere. 

n By  1971  these  programs  had  made  considerable  progress,  reducing 

the  uncertainties  in  particular  calculated  values  such  as  integrated  dose, 
to  on  the  order  of  a factor  of  two.  Also  at  this  time  the  United  States 
was  heavily  involved  in  the  design  and  deployment  of  an  anti -ballistic 
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missile  (ABM)  system.  The  Safeguard  ABM  system  consisted  of  two 
missiles,  Spartan,  an  exoatmospheric  interceptor,  and  Sprint,  which 
intercepted  reentry  vehicles  which  had  penetrated  the  atmosphere. 

Given  the  mission  of  the  latter  it  was  possible  to  conceive  of  the 
deployment  of  several  such  systems  in  close  proximity  in  time  and  space, 
with  the  resulting  possibility  of  mutual  destruction  or  as  it  is  termed, 
fratricide.  This  scenario  gave  rise  to  an  analytical  problem  in  which  the 
source,  intervening  media  and  target  were  all  known  more  or  less  exactly. 
Therefore,  the  only  uncertainties  in  analyzing  the  response  of  one  Sprint 
missile  to  the  nearby  detonation  of  another  lay  in  the  accuracy  to  which 
each  radiation  component  could  be  calculated  and  the  effect  of  each  com- 
ponent on  the  target  vehicle. 

Based  on  related  studies,  SAI  suggested  to  DNA  that  neutron  inter- 
actions within  the  target  missile  could  conceivably  produce  gamma  radia- 
tion which  was  of  the  same  order  of  intensity  as  the  gamma  environment 
impinging  on  the  missile.  Because  of  improvements  in  source  definition 
and  radiation  transport  prediction  capabilities  resulting  from  the  DNA 
programs,  it  was  both  feasible  and  desirable  to  calculate  this  (n,  missile 
y)  component.  Further,  SAI  proposed  that  the  n,  y calculations  be  made 
in  the  adjoint  mode  to  produce  a missile  radiation  response  data  base 
which  could  be  used  to  generate  internal  doses  from  any  incident  neutron 
or  gamma  ray  source.  The  data  base  would  also  generate  internal  doses 
based  on  any  response  function. 

This  effort  was  supported  by  DNA  and  led  to  results^  which  were 
of  great  interest  to  the  designers  of  the  Sprint  system,  particularly  be- 
cause for  the  first  time  precise  definitions  of  the  internal  radiation 
environment  had  been  determined.  This  new  interest  led  to  a reassess- 
ment of  the  response  of  various  key  missile  components  to  the  newly 


calculated  radiation  environments.  In  1973  at  a meeting  of  the  Basic 
Mechanisms  Panel  of  the  Transient  Effects  on  Electronic  Systems 
(TREES)  Working  Group,  several  programs  to  study  various  aspects  of 
the  response  of  missile  components  to  the  new  environments  were  pro- 
posed. Because  of  the  expense  of  these  programs  and  because  of  the 
past  shortcomings  of  radiation  transport  calculations,  the  panel  members 
were  skeptical  of  both  the  new  radiation  environments  and  the  proposed 
studies.  A particular  aspect  of  concern  was  the  adequacy  of  the  cross 
sections  and  geometric  representation  of  the  Sprint  which  were  used  in 
the  radiation  transport  calculations. 

As  a result  of  this  concern  a program  to  verify  techniques  of  cal- 
culating the  in -missile  radiation  environment  was  developed.  This  pro- 
gram consisted  of  experimental  measurements  of  a section  of  a Sprint 
missile  and  an  analytical  program  to  calculate  the  measured  data.  Four 
groups  were  involved  in  the  program.  Martin  Marietta,  the  developer 
and  manufacturer  of  the  Sprint  supplied  a detailed  description  of  the 
Sprint  missile  section  including  blueprints,  parts  lists  and  elemental 
compositions.  Inteicom  Radiation  Technology  Corporation  performed  the 
LINAC  radiation  experiments  which  were  calculated  by  Science  Applications, 
Inc.  and  Oak  Ridge  National  Laboratory. 

The  entire  program  was  formulated  as  a single  coordinated  effort 
with  each  phase  planned  jointly  by  both  experimenters  and  analysts.  In 
this  way  it  was  hoped  that  the  program  goal  could  be  accomplished  within 
the  state-of-the-art  constraints  of  both  experimental  and  analytical  cap- 
abilities. A major  problem  with  the  analysis  of  past  experiments  had 
been  in  the  definition  of  detector  response  as  a function  of  radiation  type 
and  energy  and,  perhaps  more  important,  the  communication  of  this 
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information  to  the  analyst.  This  was  seen  as  a particularly  difficult 
problem  for  this  project  since  space  limitations  dictated  the  fabrication 
of  an  entirely  new  detector  for  the  experimental  task.  As  a result  of 
these  considerations  a preliminary  experiment  was  designed  which  sub- 
jected the  detector  to  an  environment  similar  to  that  which  it  might  see 
in  the  missile  but  in  a very  simplified  geometry.  This  was  intended  to 
give  analysts  an  opportunity  to  refine  their  representation  of  the  detector 
in  the  absence  of  other  complicating  factors. 

It  was  also  felt  that,  aside  from  uncertainties  in  the  cross  section 
information  to  be  used,  the  greatest  source  of  calculational  error  would 
probably  be  inaccurate  geometric  representation  of  the  complex  components 
of  the  missile  system.  To  ensure  maximum  reliability  of  that  computational 
model  another  preliminary  experiment  was  planned.  In  this  one,  the  de- 
tectors were  placed  in  the  missile  which  was  then  exposed  to  a Cs-137 
photon  source.  This  was  intended  to  give  analysts  an  opportunity  to  check 
their  model  with  a relatively  simple  photon  only  transport  problem  utilizing 
cross  sections  more  certain  and  far  less  complex  than  is  the  case  for 
neutrons. 

Detector  placement  and  missile  orientation  in  the  final  experiment 
were  based  on  the  internal  radiation  environments  of  the  previous  SAI  cal- 
culations. Detectors  were  placed  in  the  missile  guidance  set  (MGS)  and 
the  autopilot,  and  the  missile  and  detectors  were  exposed  to  neutrons 
from  the  IRT  Linac  facility  in  a side -on  and  nearly  head-on  configuration. 

In  addition,  to  obtain  information  on  low  energy  as  well  as  high  energy 
neutrons,  the  missile  was  exposed  at  two  flight  path  stations,  17.  59 
meters  for  the  low  energy  neutron  flux  and  53.09  meters  for  the  high 
energy  neutron  flux.  Budget  considerations  dictated  that  the  final 
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analyses  of  the  experiments  be  accomplished  by  one  organization  only. 
Therefore,  SAI  analyzed  both  the  preliminary  experiments  and  the  final 
experiments  while  ORNL  calculated  only  the  preliminary  experiment. 

The  results  of  the  SAI  analysis  are  presented  in  the  following  sections. 

The  quality  of  both  the  experimental  and  calculational  results  of 
this  project  cannot  be  fully  appreciated  without  a knowledge  of  the  dramatic 
change  in  the  state-of-the-art  of  predicting  radiation  environments  occurring 
over  the  decade  from  1965  - 1975.  For  more  detailed  discussions  of  pre- 
vious experimental  attempts  to  determine  radiation  environments  a review 
of  the  literature  is  required.  For  example  during  the  1960 's  two  projects, 
Bare  Reactor  Experiment -Nevada  (BREN)  and  High  Energy  Neutron  Re- 
actions Experiment  (HENRE)  were  designed  to  obtain  dose  versus  range 
for  various  source  heights.  The  reactor  (BREN)  and  accelerator  (HENRE) 
were  placed  on  a tower  in  Nevada  and  various  measurements  were  made 
on  the  ground  and  in  the  air  by  suspending  detectors  from  helicopters. 

Also  other  experiments  were  designed  by  Gulf  General  Atomic  to  verify 
the  time  dependent  radiation  transport  in  concrete^-^  and  by  Lawrence 

Livermore  Laboratory  to  verify  cross  sections  (pulsed  sphere  experi- 
(9) 

ments  ).  Finally  it  should  be  noted  that  the  techniques  reported  here 
are  the  result  of  the  use  of  what  is  now  off-the-shelf  technology  and  not 
merely  that  of  a special  effort  designed  to  solve  a single  problem.  Thus, 
these  results  are  an  indication  of  the  degree  of  excellence  which  the 
state-of-the-art  of  radiation  transport  prediction  has  attained. 
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2.  steel-styrofoam  design  effort 


Because  of  the  geometrical  complexities  of  the  Sprint  missile,  it 
was  decided  early  in  the  program  that  a simpler  geometry  should  also  be 
measured  and  calculated.  This  "simple”  experiment  would  allow  the  ex- 
perimeters and  calculaters  to  compare  some  results  not  clouded  by  geo- 
metry or  cross  section  complexities.  In  order  for  the  simple  experiment 
to  be  a useful  comparison,  it  had  to  meet  several  criteria.  First  its  total 
and  time  dependent  response  should  be  as  much  like  the  Sprint  missile  as 
possible.  Second  its  geometry  should  be  well  defined  and  its  cross  sections 
well  known.  Third  calculation  or  measurement  of  the  simple  experiment 
should  not  present  any  additional  difficulties  not  expected  in  the  Sprint  ex- 
periment. The  fourth  and  final  criteria  was  that  the  sample  be  constructed 
at  a reasonable  cost. 

With  these  guidelines  in  mind,  SAI  was  given  the  task  of  designing 
the  simple  experiment.  Initially  a sphere  design  was  considered  as  the 
simplest  possible  geometry.  Previous  adjoint  Monte  Carlo  calculations 
by  SAI  of  the  Sprint  missile  response  function  were  folded  with  a point  silicon 
flux-to-dose  response  for  neutrons  and  gamma  rays  to  approximate  the  re- 
sponse of  the  silicon  surface  barrier  detector  inside  the  missile.  The  units 
of  this  missile  response  are  rads(Si)/(incident  neutr on/cm  ).  Thus  a neu- 
tron flux  on  the  missile  surface  is  converted  to  dose  in  silicon  which  includes 
all  radiation  transport  effects  such  as  neutron  attenuation,  inelastic  and 
capture  gamma-ray  production,  and  neutron  and  gamma-ray  scattering. 
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In  order  to  design  the  one -dimensional  spherical  model  of  the  Sprint, 
an  initial  geometry  was  selected  similar  in  size,  density  and  composition  to 
the  Sprint  Autopilot  section.  Starting  with  this  initial  configuration,  adjoint 
Monte  Carlo  calculations  with  the  MORSE  code^1*’  ^ of  the  silicon  dose  of 
several  similar  sphere  configurations  were  compared  with  the  desired 
missile  response.  The  initial  geometry  configuration  was  composed  of 
H,  C,  A1  and  Fe  so  that  the  number  of  atoms  of  each  matched  the  number 
of  similar  elements  in  the  combinatorial  geometry  model  of  the  Autopilot 
section.  The  Autopilot  section  was  found  to  have  a volume  of  1.04  x 104cm3 
excluding  voids  and  containing  the  following  number  of  atom  types: 

H 1.640  (+26) 

C 2. 04 1 (+26) 

Al  (11  <;  z <;  20)  1. 090(+26) 

Fe  (z  > 20)  . 534(+26) 

In  order  to  use  polystyrene  CjHj  as  the  source  of  hydrogen  and  carbon, 
the  excess  carbon  was  included  as  aluminum  so  that  the  assumed  number 
of  atoms  were 


CH 

1.  640  (+26) 

Al 

1.  490(+26) 

Fe 

. 534(+26) 

3 

From  the  total  volume  of  1.04(+4)  cm  , the  number  of  atoms  and  assuming 
natural  densities  for  Al  and  Fe,  radii  for  the  Fe,  Al  and  CH  and  the  CH 
density  were  calculated  as 


r(cm) 

P (gm/cm  ) 

CH 

12.04 

.484 

Fe 

12.38 

7.86 

Al 

13.54 

2.7 

13 

Adjoint  Monte  Carlo  calculations  of  the  silicon  response  to  this  and 
other  similar  spheres  were  performed  with  the  Monte  Carlo  code  MORSE. 
During  these  preliminary  calculations,  possible  estimators  for  the  final 
Sprint  calculations  were  investigated.  For  the  one  dimensional  sphere, 
the  most  efficient  estimator  was  boundary  crossing  on  the  outer  surface. 
However  a next  flight  estimator  to  an  infinite  monodirectional  beam  worked 
almost  as  well  as  the  boundary  crossing  estimator,  and  thus  was  selected 
for  the  three  dimensional  Sprint  calculations.  The  form  of  the  adjoint  next 
flight  estimator  to  the  monodirectional  beam  is 

“T  • 

score  = W j P. . f- ((?•(?')  e 3 

where  scoring  is  performed  before  each  collision  and  source  generation 
where 

W = the  particle  weight 

£?•  O ' = the  scattering  cosine  to  the  beam  direction 

P = the  downscatter  probability  but  is  not  normalized  as  it 
1 includes  effects  of  absorption  or  secondary  generation 

f. j = the  normalized  scattering  angular  distribution,  and 
r j = the  mean  free  paths  to  the  outside  surface. 

o 

This  equation  is  like  a flux -at -a -point  estimator  without  the  l/4n-R  factor. 

A forward  calculation  with  next  flight  estimation  to  the  sphere  center 
was  also  performed.  The  results  agreed  with  the  adjoint  calculations  but  be- 
cause of  the  infinite  variance  nature  of  the  estimator,  the  statistics  were 
poorer.  Thus  at  this  time  it  appeared  that  the  adjoint  calculation  with  next 


flight  estimation  to  neutron  beam  would  be  the  most  efficient  technique 
for  the  final  Sprint  experiment  calculations. 

The  first  adjoint  calculations  of  the  initial  sphere  showed  that 
thermal  neutron  buildup  in  the  CH  sphere  was  greater  than  similar  cal- 
culations of  the  Sprint  Autopilot,  most  likely  because  the  materials  are 
more  homogeneous  in  the  Sprint.  However  reducing  the  density  of  the 
CH  to  . 34  gm/cc  or  .01575  atoms/barn-cm  each  of  H and  C and  changing 
the  CH  radius  to  12. 4 cm  gave  a similar  number  of  scatterings  as  the 
Sprint  Autopilot  calculations.  Even  so  with  this  new  geometry,  buildup 
caused  enough  statistical  fluctuation  that  calculations  of  less  than  one 
minute  of  the  CDC  7600  were  not  adequate  to  optimize  the  design  of  the 
outer  metal  shell. 

In  order  to  optimize  the  design  calculations  of  the  outer  material, 
a simple  correlated  sampling  version  MORSE  was  developed  so  that  several 
thicknesses  of  the  outer  material  could  be  investigated  with  the  same  random 
walks.  Several  particle  weights  were  kept  for  each  history,  each  weight 
corresponding  to  a different  thickness  of  the  outer  material,  and  the  outer 
material  was  modeled  as  several  concentric  spheres.  Whenever  a history 
crossed  one  of  these  spheres  in  an  outward  direction,  the  corresponding 
weight  was  scored  in  a leakage  estimator  and  then  that  weight  was  set  to 
zero  so  that  if  the  same  history  later  recrossed  the  same  surface  it  would 
score  zero.  This  simple  correlated  sampling  trick  allowed  several  outer- 
layer  thicknesses  to  be  investigated  on  a single  calculation,  and  because 
they  used  the  same  random  histories  the  results  were  correlated.  This 
means  that  differences  between  thicknesses  are  not  related  to  the  absolute 
uncertainties  in  the  Monte  Carlo  results  and  thus  are  a better  measure  of 
the  effect  of  the  additional  thickness  than  would  be  obtained  by  several  in- 
dependent calculations. 


Outer  materials  that  were  investigated  were  Al,  Fe,  and  Cu.  It 
was  found  that  all  three  compared  adequately  with  the  Sprint  response  for 
high  energy  neutrons,  but  at  energies  below  500  keV  both  Al  and  Fe  were 
much  lower  than  the  desired  response.  Table  1 compares  the  responses 
of  Al,  Cu,  and  Fe  at  three  radii  to  the  desired  response.  Apparently  the 
missile  response  at  low  energies  is  high  because  of  the  large  eV  capture 
resonance  in  tantalum  which  is  located  close  to  the  calculation  point  in  the 
missile.  Copper  has  a similar  resonance  at  229  eV,  although  not  as  strong 
as  tantalum,  and  thus  the  2. 1 cm  thickness  of  copper  gives  the  response 
closest  to  the  desired  response.  However  at  an  early  meeting  of  experi- 
menters and  calculators,  it  was  decided  to  construct  the  experiment  with 
iron  instead  of  copper  because  the  energy  spectra  of  the  epithermal  capture 
gamma  rays  are  better  known  for  iron.  Also  copper  activation  posed  pos- 
sible experimental  difficulties.  Dr.  Larry  Harris  of  IRT  also  requested 
that  the  sample  be  rectangular  or  cylindrical  rather  than  spherical  to 
facilitate  fabrication. 

At  this  point  several  steel-styrofoam  configurations  were  calculated 
to  find  an  easily  fabricated  design  that  would  not  have  any  calculation  or 
experimental  difficulties.  Each  of  these  calculations  had  approximately  the 
same  volume,  iron  mass  and  CH  mass,  the  main  differences  being  geometry. 
Rather  than  documenting  each  of  these  calculations  separately,  they  are 
listed  in  Table  2 along  with  the  reason  for  the  calculation  and  a brief  con- 
clusion. Because  these  geometries  did  not  have  spherical  symmetry,  the 
leakage  estimator  used  in  the  previous  correlated  sampling  sphere  calcu- 
lations was  no  longer  valid.  Each  of  these  calculations  used  the  adjoint 
next  flight  estimator  to  the  parallel  beam  which  was  briefly  discussed 
earlier  in  this  chapter.  When  a next  flight  estimator  is  used  to  compute 
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TABLE  2.  STEEL-STYROFOAM  DESIGN  CALCULATIONS 
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TABLE  2 


INCIDENT  BEAM  INCIDENT  BEAM 


Figure  1.  Wafer  Design  for  Steel  Styrofoam  Sample 
(Calculation  4 on  Table  2) 


Figure  2.  Double  Wafer  Design  for  Steel  Styrofoam 
Sample  (Calculation  5 on  Table  2) 


the  flux  at  a point  in  a forward  calculation,  an  infinite  variance  condition 

2 

often  arises  because  of  the  (1/4 ffr  ) factor  in  the  score  when  particles 
scatter  near  the  point,  causing  large  variances  that  are  not  readily  re- 
duced by  calculating  more  histories.  This  problem  does  not  arise  with 

the  adjoint  next  flight  estimator  to  a beam  because  the  estimation  to  a 

2 

beam  has  no  (l/4ffr  ) factor  and  no  scatterings  near  an  estimation  point. 

From  Table  2 the  important  features  governing  the  design  of  the 
cylinder  experiment  can  be  seen.  From  1,  2,  and  6 compared  with  4 and 
5 it  is  important  that  the  steel  entirely  surround  the  styrofoam  in  order  to 
maximize  leakage.  Low  energy  neutrons  are  thermalizing  in  the  CH  and 
then  being  captured  in  the  Fe  with  the  detector  seeing  the  capture  gamma 
rays.  To  maximize  the  dose  from  low  energy  neutrons,  they  must  not  be 
able  to  leak  out  of  the  CH  without  entering  the  Fe.  This  is  the  reason  that 
the  waffer  designs  of  4 and  5 (Figs.  1 and  2)  give  such  a low  response. 

The  box  (calculation  2)  gave  good  results  and  would  be  easiest  to  fabricate, 
but  calculation  3 showed  that  the  beam  fringes  acting  on  the  corners  would 
have  significant  effects.  Because  the  beam  edges  cannot  be  exactly  deter- 
mined, the  cylinder  design  was  adopted  to  fit  the  entire  sample  into  the  beam. 
Even  so,  calculation  7 showed  that  . 8 cm  difference  in  the  beam  diameter 
could  have  50%  effects  on  the  response  to  neutrons  from  . 5 MeV  to  3 keV. 

For  this  reason  care  was  taken  to  assure  that  the  entire  steel  cylinder  fit 
safely  into  the  beam.  Calculation  8 showed  no  effects  of  the  incident  beam 
striking  the  cylinder  walls  at  a small  angle  and  calculation  9 showed  that 
the  steel  thickness  chosen  by  correlated  sampling  of  the  sphere  sample 
was  still  correct  for  the  cylinder. 

The  dimensions  of  the  cylinder  used  in  the  above  calculation  were 
constrained  so  that  (1)  the  total  volume  would  be  the  same  as  the  previous 
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correlated  sphere  calculation,  (2)  the  cylinder  diameter  was  20.  32  cm, 
and  (3)  the  iron  thickness  was  1.  6 cm.  These  constraints  gave  the  cylinder 
a length  of  35.  44  cm.  For  final  construction  at  Intelcom  the  length  was 
reduced  to  the  diameter  to  make  a more  regular  sample.  Figure  3 shows 
the  final  dimensions  of  the  steel  styrofoam  cylinder  constructed  by 
Intelcom  and  used  in  their  experiment. The  cylinder  and  end  plates 
were  made  of  slightly  different  steels. 


(Dimensions  in  cm) 


Figure  3.  Sectional  View  of  the  Steel  Styrofoam 
Sample  (Furnished  by  IRT) 
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3.  RADIATION  TRANSPORT  MODEL 


The  heart  of  any  three  dimensional  Monte  Caclo  calculation  is  the 
radiation  transport  model  which  consists  of  (1)  a geometry  model  depicting 
the  missile  components  to  the  desired  accuracy,  (2)  the  elemental  den- 
sities of  each  region  modeled  in  the  geometry,  and  (3)  the  physical  cross 
sections  for  the  elements  used  in  the  model.  The  first  section  of  this 
chapter  presents  our  sources  of  cross  section  data,  and  describes  our 
efforts  necessary  to  process  this  data  into  a multigroup  library  compatible 
with  the  MORSE  multigroup  Monte  Carlo  code.  The  next  two  sections 
document  both  the  steel-styrofoam  model  and  the  Sprint  missile  model 
respectively.  Data  for  the  steel -styrofoam  sample  was  supplied  by  LRT 
and  the  Sprint  was  modeled  from  data  supplied  by  Martin -Marietta  Corp. 
The  final  section  describes  the  neutron  source  and  the  dosimeter  response 
which  were  supplied  by  IRT. 

3.  1 CROSS  SECTIONS 

Coupled  fine  group  neutron -gamma  cross  section  sets  were  pro- 
duced by  processing  ENDF/B  and  DNA  cross  section  libraries  with  the 
(13) 

AMPXV  system  of  codes  which  is  briefly  described  below.  The  fine 
group  structure  that  was  selected  is  described  in  Tables  3 and  4.  These 
in  turn  were  collapsed  into  broad  groups  using  the  code  COLIXS  and 
weighting  functions  produced  by  XSDRNPM--a  module  of  AMPX.  The 
broad  group  structure  is  presented  in  Table  5 along  with  the  dosimeter 
response.  The  elements  listed  in  Table  6 have  been  processed  and  are 
available  in  both  the  fine  and  the  broad  group  structure. 
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TABLE  3.  NEUTRON  FINE  GROUP  STRUCTURE 


Group Energy  Range Lethargy  Range 
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1 . 9 .3  0 0 7 4 0 0 

- 

2 , 4 7 9 0 p 404 

6.000 

6.250 

63 

1 ,S050t 400 

- 

1 .930074  04 

6.251 

6.500 

b'l 

1 , 1 7 1 0 F 4 o O 

- 

1 .50307+04 

6 . 50  n 

- 

6.750 

65 

9, 1 1 40F  40  3 

- 

1 . 17104  *oa 

6.750 

- 

7,000 

66 

7, 1 0 ? 0 E 4 0 5 

- 

9.11907403 

7.000 

- 

7.26  0 

67 

5.5310'  403 

- 

7.  1 P ? o 7 4 o 3 

7 . ? 5 0 

- 

7.500 

60 

0, 5 0 7 o F 403 

- 

5.53107  4 13 

7.504 

- 

7.750 

69 

5.5550F403 

- 

4. 30704  4 03 

7.750 

- 

6,000 

70 

2.61307 403 

- 

3.35507*03 

6.000 

- 

6.250 

7 1 

2. 0 350F  403 

- 

2.61 507403 

6.250 

- 

6.500 

7? 

1 .SS50F40  1 

- 

2.0  35OP  «03 

6.500 

- 

6.760 

73 

1 , ? mop  f 03 

- 

1,56507403 

6,750 

- 

9.000 

70 

9,61  10*  *02 

- 

1 .234074  03 

9,000 

- 

9.260 

75 

7 , 46607  4 0? 

- 

9,61  1 04  4 0? 

9.250 

9 . 5 0 0 

76 

5.6290540? 

- 

7.44507+02 

0,500 

- 

9.750 

7’ 

3 . 5 •,  0 0 F 4 r ? 

• 

5.62907+02 

9,750 

- 

10.000 

7 6 

5 , 536 or *0? 

- 

4,54007 + n ? 

10.000 

10.250 

/9 

?.  75404  40? 

- 

5,5360740? 

10.250 

- 

10.500 

«o 

?,  1450F40? 

- 

2.  76407  4 0? 

10.500 

- 

10.750 

61 

1 . 6 7 0 0 7 40? 

- 

2,  1 4504  4 0? 

10.750 

- 

11.000 

62 

1.501  OF  4 0? 

- 

1 .67007  4 0? 

1 1 .OOO 

- 

1 1 .250 

65 

1,01304 40? 

- 

1.5010240? 

1 1 .250 

1 1 .500 

64 

7,5f>°C7  40  1 

- 

1.0150440? 

1 1,500 

1 1 . 750 

65 

6 , 1 4 4 0 F + 0 I 

- 

7.6490740' 

1 1 .750 

12.000 

66 

4 . 7650E  + 0 1 

- 

6. 1 4404  *0  1 

12.000 

1 2.250 

67 

5,72707461 

- 

4,76507401 

12.250 

- 

1 2.500 

66 

2, 9020^40 | 

- 

3,72/07401 

12.500 

12.760 

6') 

?,?600K401 

- 

2,90204  *01 

12.750 

- 

1 3 . 0 n 0 

90 

1 ,7  60  04  40  1 

- 

2.26002*01 

13,000 

1 3.260 

91 

1 . 57  1 Of  401 

- 

1 .76007*01 

13.250 

- 

1 3.500 

9? 

1 ,06807  401 

- 

1 .37101  *01 

13.500 

- 

1 3.  /50 

95 

6, 5150F4O0 

• 

1 .06604  *0  1 

13.756 

- 

14.000 

90 

6,4  76  04  4 0 0 

- 

6.  515074  on 

14.000 

- 

1 n,?60 

95 

5 . 0 4 5 0 fc  4 0 0 

- 

6.47607  400 

14.250 

- 

1 0.500 

96 

5, 9260 F 4. 10 

- 

5,00307400 

14.500 

- 

1 n . 760 

97 

5,05904  400 

- 

5.9?(«<'7*00 

14.750 

15.000 

96 

2.36204 400 

- 

5.05907  4 00 

15.000 

- 

15.260 

99 

1 ,65507400 

- 

2 . 36202 « oo 

15.250 

- 

15.500 

1 on 

J , 40507400 

- 

1 .05504 400 

15.600 

“ 

15, 750 

101 

1 .,12507400 

- 

1 .04507400 

15.750 

- 

16.000 

102 

6 .‘7640 7 -o  1 

• 

1 . 1 2507  400 

16.000 

- 

16.250 

105 

6.62604 -o 1 

- 

6.76401 -01 

16.250 

16.600 

ion 

5,51607-01 

- 

6. 62607-01 

16.500 

- 

16,760 

105 

0,14007-01 

- 

5.31607-01 

16.750 

- 

17.000 

1 06 

1,00007-05 

- 

4.14007-01 

17.000 

- 

27,631 
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TABLE  4.  GAMMA-RAY  FINE  GROUP  STRUCTURE 


Group  Enorpy  Range 


1 

1.000  E 

+ 

07 

- 

1.200  E 

+ 07 

2 

9. 000  E 

+ 

06 

- 

1.000  E 

+ 07 

3 

8.000  E 

+ 

06 

- 

9.000  E 

+ 06 

4 

7. 500  E 

+ 

06 

- 

8.000  E 

+ 06 

5 

7. 000  E 

+ 

06 

- 

7.500  E 

4 06 

6 

6.500  E 

+ 

06 

- 

7.000  E 

+ 06 

7 

6.000  E 

+ 

06 

- 

6.500  E 

4-  0Q 

8 

5.500  E 

+ 

06 

- 

6.000  E 

4-  05 

9 

5.000  E 

+ 

06 

- 

5. 500  E 

4-  06 

10 

4.500  E 

+ 

06 

- 

5.000  E 

4-  06 

11 

4.000  E 

+ 

06 

- 

4. 500  E 

+ 06 

12 

3.500  E 

+ 

06 

- 

4.000  E 

+ 06 

13 

3.000  E 

+ 

06 

- 

3. 500  E 

4-  06 

14 

2.500  E 

+ 

06 

- 

3.000  E 

4 06 

15 

2.000  E 

+ 

06 

- 

2.500  E 

4-  06 

16 

1.  660  E 

+ 

06 

- 

2.000  E 

4-  06 

17 

1.500  E 

+ 

06 

- 

1.  660  E 

4-  06 

18 

1.330  E 

+ 

06 

- 

1.500  E 

4-  06 

19 

1.000  E 

+ 

06 

- 

1.330  E 

4-  06 

20 

8.000  E 

+ 

05 

- 

1.000  E 

+ 06 

21 

7.000  E 

+ 

05 

- 

8.030  E 

+ 05 

22 

6.000  E 

+ 

05 

- 

7.000  E 

4 05 

23 

5.000  E 

+ 

05 

- 

6.000  E 

+ 05 

24 

4.000  E 

+ 

05 

- 

5.000  E 

4 05 

25 

3.000  E 

+ 

05 

- 

4.000  E 

+ 05 

26 

2.000  E 

+ 

05 

- 

3.000  E 

4-  05 

27 

1.000  E 

+ 

05 

- 

2.000  E 

+ 05 

28 

7. 500  E 

+ 

04 

- 

1.000  E 

4-  05 

29 

5.000  E 

+ 

04 

- 

7.  500  E 

4-  0 4 

30 

3.000  E 

+ 

04 

- 

5.000  E 

4-  04 

31 

1.000  E 

4* 

04 

- 

3.000  E 

+ 04 
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TABLE  5.  BROAD  GROUP  COUPLED  NEUTRON  AND 

GAMMA  RAY  ENERGY  BOUNDS  AND  SILICON 
DOSIMETER  RESPONSE 


Group 

Upper  Energy  Bound 
(eV) 

Silicon  Dosimeter 
Response 

rad(Si)/(n  or  y/cnri  ) 

1 

1.50  (+7) 

1.  15(-9) 

2 

1.  35(+7) 

1.  02(-9) 

3 

1.  U(+7) 

8.  65(-10) 

4 

9. 05(+6) 

6.  1 7( - 10) 

5 

7. 00(+6) 

2.  77 (-10) 

6 

5.  45(+6) 

1.  36(-10) 

7 

4. 10  (+6) 

7.  40 (-11) 

8 

3.01(+6) 

4.  70 (-11) 

9 

2. 41(+6) 

3.  90(-ll) 

10 

2.  30 (+6) 

2.  71 (-1 1) 

11 

1.  85(+6) 

1.  23(-ll) 

12 

1.  ll(+6) 

0. 

13 

5. 50(+5) 

0. 

14 

1.  60(+5) 

0. 

15 

1.  ll(+5) 

0. 

16 

2.  50 (+4) 

0. 

17 

1.  20  (+4) 

0. 

18 

3. 40(+3) 

0. 

19 

5.  90(+2) 

0. 

20 

1.01  (+2) 

0. 

21 

1.  80(+l) 

0. 

22 

3.  10 

0. 

23 

4. 14(-1) 

0. 

24 

1.  20(+7) 

2. 40 (-9) 

25 

8. 00(+6) 

1.  90 (-9) 

26 

6.  50(+6) 

1.  60 (-9) 

27 

5. 00(+6) 

1.  30(-9) 

28 

3.  50 (+6) 

9.  70 (-9) 

29 

2.  50 (+6) 

6.  90(-10) 

30 

1.50  (+6) 

4.  60(-10) 

31 

1.00  (+6) 

4.  60(-10) 

32 

7. 00  (+5) 

3.  40(-10) 

33 

4. 00  (+5) 

2.  20(-10) 

34 

2. 00(+5) 

8.  70(-ll) 

35 

7.  50  (+4) 

0. 

36 

5. 00(+4) 

0. 
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TABLE  6.  MATERIALS  PROCESSED  BY  AMPX  FOR  USE 
IN  THE  SPRINT  CALCULATIONS 


Material 

Source  of  Cross  Section  (x/s)  Data 

H-l 

DNA  4184  Mod  2 - thermal  x/s  adjusted  for  bound  H 

H-3 

DNA  41C9  Mod  1 - no  gamma  production 

He -4 

DNA  4504  Mod  0 - no  gamma  production 

Be -9 

DNA  4154  Mod  2 

B-10 

ENDF/B  1155  - gamma  production  x/s  added  for  n,  cry 
only 

C 

DNA  4274  Mod  0 

N 

DNA  4133  Mod  4 

O 

DNA  4134  Mod  2 

F 

ENDF/B  919  for  neutron  x/s,  DNA  4509  for  gamma 
production 

Na 

DNA  415G  Mod  1 

Mg 

DNA  4512  Mod  1 

Ai 

DNA  4135  Mod  3 

Si 

DNA  4151  Mod  3 

Cl 

ENDF/B  1149 

Ca 

DNA  4152  Mod  3 

Ti 

ENDF/B  1286 

Cr 

ENDF/B  1191 

Fe 

DNA  4180  Mod  2 

Ni 

ENDF/B  1190 

Cu 

DNA  4529  Mod  1 

Zn 

DNA  4529  Mod  1 (Cu)  for  neutron  x/s,  DNA  4530  (Zn)  for 
gamma  production 

Zircaloy 

ENDF/B  1284  for  neutron  x/s,  ENDF/B  1189  (Nb)  for 
gamma  production 

29 


T 


TABLE  6.  (Cont'd.) 


< 


Material 

Source  of  Cross  Section  (x/s)  Data 

Nb 

ENDF/B  1189 

Ag-107 

ENDF/B  1138  for  neutron  x/s,  DNA  4547  for 
production 

gamma 

Ag-109 

ENDF/B  1139  for  neutron  x/s,  DNA  4547  for 
production 

gamma 

Cd-NAT 

ENDF/B  1281  for  neutron  x/s,  ENDF/B  1287 
gamma  production 

(Mo)  for 

Sn 

ENDF/B  1138  (Ag-107)  for  neutron  x/s,  DNA 
for  gamma  production 

4550  (Sn) 

Gd-NAT 

ENDF/B  1030  for  neutron  x/s,  ENDF/B  1291 
gamma  production 

(Eu)  for 

Ta-181 

DNA  4179  Mod  3 

W-182 

DNA  4582  Mod  2 

W-183 

DNA  4583  Mod  3 

W-184 

DNA  4584  Mod  3 

W-186 

DNA  4586  Mod  2 

Au-197 

DNA  4283  Mod  0 

Pb 

DNA  4136  Mod  5 

U-235 

DNA  4188  Mod  1 

U-238 

DNA  4187  Mod  1 

Pu-239 

DNA  4539  Mod  0 

Pu-240 

DNA  4540  Mod  0 
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Cross  sections  for  the  Cs  137  calculations  were  processed  with 
the  SMUG  module  of  the  AMPX  code  for  the  elements  shown  in  Table  7 
into  a fine  group  structure  and  collapsed  to  the  broad  groups  shown  in 
Table  8.  The  source  of  the  x-ray  cross  sections  is  the  Radiation  Shield- 
ing Information  Center  DLC-7  library. 

AMPX  is  a modular  system  for  producing  coupled  multigroup 
neutron -gamma  cross  section  sets.  Basic  cross  section  data  for  AMPX 
are  obtained  from  ENDF/B  like  libraries.  Most  commonly  used  opera- 
tions required  to  generate  and  collapse  multigroup  cross  section  sets  are 
provided  in  the  system.  AMPX  was  assembled  at  ORNL  by  N.  M.  Greene, 
et  al.  for  use  on  the  IBM  360  series  machines.  SAI  has  converted  the 
code  for  use  on  CDC  7600  and  6600  machines. 

The  basic  processes  provided  will  (1)  generate  multigroup  neutron 
cross-sections;  (2)  generate  multigroup  gamma  cross  sections;  (3)  gene- 
rate gamma  yields  i ir  gamma-producing  neutron  interactions;  (4)  combine 
neutron  cross  sections,  gamma  cross  sections,  and  gamma  yields  into 
final  "coupled  sets";  (5)  perform  one -dimensional  discrete  ordinates  trans- 
port or  diffusion  theory  calculations  for  neutrons  and  gammas  and,  an 
option,  collapse  the  cross  sections  to  a broad-group  structure,  using  the 
one  -dimensional  results  as  weighting  functions;  and  (6)  output  multigroup 
cross  sections  in  convenient  formats  for  other  codes. 

AMPX  includes  the  following  modules. 

1.  DRIVER  - The  AMPX  module  manager. 

2.  XLACS  - Produces  weighted  multigroup  neutron  cross  sections 
from  ENDF/B. 

3.  LAPHNGAS  - Generates  multigroup  secondary  gamma-ray 
production  cross  sections. 
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TABLE  7.  ELEMENTS  PROCESSED  INTO  THE  X-RAY  GROUP 
STRUCTURE  FOR  THE  CS  137  CALCULATION 


Atomic  Number 

Element 

1 

H 

2 

He 

3 

Li 

4 

Be 

5 

B 

6 

C 

7 

N 

8 

O 

9 

F 

11 

Na 

12 

Mg 

13 

A1 

14 

Si 

17 

Cl 

18 

Ar 

20 

Ca 

22 

Ti 

24 

Cr 

26 

Fe 

28 

Ni 

29 

Cu 

30 

Zn 

42 

Mo 

47 

Ag 

50 

Sn 

56 

Ba 

73 

Ta 

74 

W 

79 

Au 

82 

Pb 

92 

U 

94 

Pu 
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TABLE  8.  X-RAY  GROUP  STRUCTURE  FOR  THE  CS  137 

CALCULATIONS  AND  THE  ASSUMED  DOSIMETER 
RESPONSE 


Group 

Upper  Energy 
(keV) 

Silicon  Dose  ? 
rad(Si)/(x/cm  ) 

1 

662. 

2.95  (-10) 

2 

600. 

2.69  (-10) 

3 

500. 

2.20  (-10) 

4 

425. 

1.83  (-10) 

5 

350. 

1.47  (-10) 

6 

275. 

1.16  (-10) 

7 

220. 

8.93  (-11) 

8 

160. 

7.22  (-11) 

9 

120. 

7.22  (-11) 

10 

90. 

9.8  (-11) 

11 

65. 

0 

12 

45. 

0 

10. 

4.  CHOX  - Performs  cross  section  interface  management. 

5.  SMUG  - Generates  multigroup  photon  cross  sections. 

6.  N1TAWL  - Provides  for  resonance  self -shielding  and  work- 
ing library  production. 

7.  XSDRNPM  - Provides  one -dimensional  Sn  capability  for 
spatial  weighting. 


3.  2 STEEL  STYROFOAM  MODEL 

The  steel  styrofoam  experiment  was  modeled  by  six  regions,  three 
of  which  modeled  the  dosimeter  and  were  used  in  the  Sprint  model  as  well. 
Figure  4 shows  the  six  regions  of  the  steel  styrofoam  experiment  which 
are  as  follows: 

1.  surface  barrier  detector 

2.  surrounding  silicon  and  Al  can 

3.  attaching  wires  and  connectors 

4.  styrofoam 

5.  front  and  back  steel  plates 

6.  steel  cylinder 

Table  9 lists  the  computer  input  necessary  to  model  the  steel  styrofoam 
experiment  with  the  combinatorial  geometry.  Combinatorial  geometry  is 
a ray  tracing  modual  of  MORSE  that  includes  all  space  as  intersections  and 
unions  of  simple  body  types  such  as  spheres,  cylinders,  rectangles,  cones 
and  others.  The  first  section  of  the  table  defines  the  necessary  bodies, 
and  the  following  cards  describe  the  intersections  and  unions  of  the  bodies. 
Table  10  lists  the  elemental  compositions  of  each  of  the  six  regions  used 
in  our  calculation. 
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ALL  DIMENSIONS  IN  CM 


Figure  4.  Six  Region  Combinatorial  Geometry  Model 
of  the  Steel  Styrofoam  Sample 
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3.  3 SPRINT  MISSILE  MODEL 


The  Sprint  missile  model  is  considerably  more  detailed  using  127 
bodies  to  generate  69  zones,  66  of  which  have  different  elemental  com- 
positions. The  model  is  the  result  of  an  evolution  which  began  in  1969 
with  an  eight  region  model  generated  by  Martin  Marietta  and  through 
various  programs  at  Martin  and  SAI  has  progressed  in  detail  to  the  exist- 
ing model.  Of  the  sixty -nine  regions,  46  are  used  to  model  the  Autopilot 
(A/P),  11  model  the  Missile  Guidance  Set  (MGS),  5 model  the  ballast  con- 
tainer and  antenna,  and  the  remainder  model  the  dosimeter  and  the  skin. 
Final  data  for  the  model  was  supplied  by  Ted  Bohm  of  Martin  Marietta  in 

(14) 

the  form  of  a letter  report v and  several  blueprints.  The  letter  report 

included  their  model  to  date,  plus  component  lists  and  compositions  of  all 
regions  in  the  vicinity  of  the  detectors.  Included  with  the  component  lists 
were  the  elemental  compositions  of  each  of  the  components. 

In  order  to  place  the  dosimeter  into  the  Autopilot  and  the  MGS,  a 
3.  41  cm  diameter  hole  was  drilled  into  each,  from  the  bottom  of  the  Auto- 
pilot and  to  the  side  of  the  MGS.  The  exact  position  of  these  holes  showing 
the  material  removed  is  detailed  in  the  Martin  letter  report.  Upon  investi- 
gating these  holes,  it  was  determined  that  the  material  removed  was  re- 
presentative of  the  homogenization,  containing  no  significant  components 
with  different  atomic  densities.  Therefore  the  Martin  Marietta  elemental 
densities  were  retained  for  the  regions  surrounding  the  detector  holes. 

The  Sprint  missile  model  used  in  our  calculations  is  documented 
in  the  next  three  tables.  Table  11  lists  the  zones  with  a brief  description, 
the  volume  and  the  mass.  Table  12  presents  the  combinatorial  geometry 
computer  input  of  the  model  and  Table  13  presents  the  elemental  com- 
position of  each  region. 
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TABLE  11.  SPRINT  MISSILE  MODEL 


Zone 

Missile 

Section 

3 

Contents  of  Region  Volume  (cm  ) 

Weight  (gm) 

1 

A/P 

Squibs 

242. 700 

354.39 

2 

A/P 

Squibs 

93. 547 

71.  98 

3 

A/P 

Accel.  Summing  Ampl. 
& Rate  Summing  Ampl. 

63.732 

70.77 

4 

A/P 

Outer  Loop  Shaping  Ampl. 

82.  744 

94.  14 

5 

A/P 

Rate  Feedback  Modulator 

190.  628 

207.07 

6 

A/P 

Inner  Loop  Shaping  Ampl. 

89.833 

96.31 

7 

A/P 

Inner  Loop  Shaping  Ampl. 

81.  940 

95.  68 

8 

A/P 

Outer  Loop  Shaping  Ampl. 
& Resistor 

85.707 

93.48 

9 

A/P 

Structural  Filter 

218.  698 

238.54 

10 

A/P 

Roll  Filter  Amplifier 

98.  949 

132. 10 

11 

A/P 

Accel.  Summing  Ampl. 
& Rate  Summing  Ampl. 
& Resistors 

78. 708 

80.99 

12 

A/P 

Dynamic  Pressure  Gain 
Change^ 

52.  161 

48.21 

13 

A/P 

T3  Timer 

59.344 

62.91 

14 

A/P 

Pulse  Width  Modulator 
and  Resistor 

69.  547 

51.41 

15 

A/P 

Area  Gain  Changer  & 
Roll  Tachometer  Ampl. 

74. Ill 

82. 16 

16 

A/P 

Connector 

51.576 

127.00 

17 

A/P 

Frame  (Areas  not  covered 
by  Regions  1-16) 

524. 772 

891.37 
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TABLE  11.  (Cont'd.) 


Zone 

Missile 

Section 

3 

Contents  of  Region  Volume  (cm  ) 

Weight  (gm) 

18 

A/P 

2nd  Stage  Servo  Ampl. 

37.605 

58.40 

19 

A/P 

Adaptive  Sig.  Inverter, 

2nd  Stage  Act.  Sum.  Ampl. 
2nd  Stage  Servo  Ampl. , 
Servo  Pwr  Ampl. , 1st 
Stage  Act.  Sum.  Ampl., 

1st  Stage  Servo  Ampl. 

9 

219.  508 

305.77 

20 

A/P 

1st  Stage  Act.  Sum.  Ampl. 
& 1st  Stage  Servo  Ampl. 

& Jumper 

55.005 

76. 12 

21 

A/P 

Comparator 

25.439 

33.90 

22 

A/P 

Adaptive  Gain  Change 
Amplifier 

38.343 

48.86 

23 

A/P 

Roll  Modulator,  Low 
Pass  Filter,  Adaptive 
Sig.  Inverter,  2nd  Stage 
Servo  Ampl. 

138. 650 

182.81 

24 

A/P 

High -Pass,  Low-Pass 
Filter,  Comparator,  1st 
Stage  Servo  Ampl. 

169.891 

211.33 

25 

A/P 

Adaptive  Integrator  & 
Jumper 

39.690 

44.  18 

26 

A/P 

Adaptive  Gain  Change 
Amplifier 

36.561 

53.48 

27 

A/P 

Signal  Control  Mod. , 2nd 
Stage  Act.  Sum.  Ampl., 
2nd  Stage  Servo  Ampl. , 
Servo  Pwr  Ampl. , 1st 
Stage  Act.  Sum.  Ampl. , 
1st  Stage  Servo  Ampl. 

166.  679 

249.  37 

40 
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TABLE  11.  (Cont’d.) 


Zone 

Missile 

Section 

Contents  of  Region 

3 

Volume  (cm  ) 

Weight  (gm) 

28 

A/P 

Adaptive  Integrator 

32. 624 

49.  36 

29 

A/P 

2nd  Stage  Act.  Sum. 
Ampl. , 2nd  Stage  Servo 
Ampl. , Servo  Pwr  Ampl, 

102.253 

140.25 

30 

A/P 

Low  Pass  Filter,  Servo 
Pwr  Ampl. 

101.571 

120.96 

31 

A/P 

Convertor 

110.451 

128.  65 

32 

A/P 

Frame  (Areas  not 
covered  by  Regions 
18-32) 

439.55 

937. 10 

33 

A/P 

Electrical  Linear 
Accelerometer 

134.06 

113.24 

34 

A/P 

Rate  Integrating 
Gyroscope 

212.14 

343.07 

35 

A/P 

Electrical  Linear 
Accelerometer 

237. 16 

194.79 

36 

A/P 

Rate  Integrating 
Gyroscope 

303.03 

324.45 

37 

A/P 

Primary  Wet  Battery 

1728. 1 

2920.  94 

38 

A/P 

Rate  Control  Elec. 
Ampl. 

207.221 

170.71 

39 

A/P 

Elec.  Control  Ampl. 
Platform 

524. 802 

765.85 

40 

A/P 

Static  Pwr  Inverter 
Subassay. 

468.359 

366.87 

41 

A/P 

Power  Ampl.  Inverter 
Subassay. 

781.724 

608.58 

42 

A/P 

Single  Axis  Platform 

892.  65 

2300.  64 

43 

A/P 

Electrical  Control  Rate 
Ampl. 

231.769 

317.  78 

TABLE  11.  (Cont'd.) 


Zone 

Missile 

Section 

Contents  of  Region 

3 

Volume  (cm  ) 

Weight  (gm) 

44 

A/P 

Capacitor  Assembly 

44.647 

8.70 

45 

A/P 

DC  Voltage  Regulator 

307.269 

1037.56 

46 

A/P 

Interconnecting  Box 

1300. 937 

1650.38 

47 

A/P 

MGS 

Dosimeter  (2) 

20. 548 

52.91 

48 

A/P 

MGS 

Dosimeter  Wires  (2) 

76.04 

11.68 

49 

Nose 

Adjustible  Ballast 
Holder 

456.  14 

907.2 

50 

MGS 

Signal  Data  Converter 
& Pulse  Decoder 

3367.0 

4326. 

51 

MGS 

Radio  Frequency 
Detector 

1642.0 

1788. 

52 

MGS 

Electrical  Frequency 
Converter 

4305.0 

4018. 

53 

MGS 

Transmitter 

2423.0 

2581. 

54 

MGS 

Transmitter 

1452.0 

799. 

55 

MGS 

Electrical  Frequency 
Converter 

117.7 

220. 

56 

MGS 

Power  Supply  & P.  S. 
Sub -Assembly 

598.8 

539. 

57 

MGS 

Power  Supply 

1825.0 

2699. 

58 

MGS 

"S"  Band  Cable  Trough 

444.5 

101. 

59 

MGS 

"L"  Band  Cable  Trough 

444.5 

156. 

60 

MGS 

Cover  & Housing 

1334.0 

4147. 

42 


TABLE  11.  (Cont’d. ) 


Zone 

Missile 

Section 

Contents  of  Region 

3 

Volume  (cm  ) 

Weight  (gm) 

61 

Skin 

Aluminum 

5480. 

5035. 

62 

Skin 

Phenolic 

4455. 

6577. 

63 

Nose 

Ballast  Container 

268. 1 

2172.2 

64 

Nose 

Ballast 

370.  36 

3079.9 

65 

Nose 

Antenna 

443.6 

1367.2 

66 

Nose 

Nose  Structure 

3810.79 

7307.5 

\ , 
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TABLE  12.  (Cont’d.) 
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Our  Sprint  model  differs  in  several  small  ways  from  the  Martin 
Marietta  model.  First  our  model  includes  the  ballast,  antennas  and  nose 
structure  which  are  unique  to  this  experiment.  Further,  several  of  our 
zones  have  slightly  different  volumes,  sometimes  more  realistic  and  some- 
times less  realistic  than  the  Martin  model  depending  upon  the  proximity  of 
the  detectors  and  other  geometry  knowledge  gained  from  other  SAI  Sprint 
programs.  However  in  all  cases  we  maintained  the  masses  documented  by 
Martin  causing  some  of  our  densities  to  be  slightly  different. 

Of  interest  in  the  conclusion  is  the  fact  that  although  46  regions 
were  used  to  model  the  autopilot,  only  eleven  were  used  for  the  MGS. 

A difference  in  the  experimental  calculation  comparisons  might  be  an 
indication  of  the  degree  of  modeling  accuracy  required  for  future  internal 
dose  calculations,  especially  if  the  MGS  results  agree  as  well  with  the 
experiment  as  the  A/P  results. 

3. 4 NEUTRON  SOURCES  AND  DOSIMETER  RESPONSE 

This  section  documents  the  neutron  sources  and  the  dosimeter 
response  which  were  supplied  by  IRT  and  used  in  our  calculations.  In 
order  to  gain  good  experimental  statistics  without  swamping  the  detector's 
counting  ability,  it  was  necessary  to  perform  the  experiment  at  two  dif- 
ferent flight  path  positions.  With  the  experiment  at  53. 095  meters  in- 
cident neutrons  from  20  MeV  to  . 18  MeV  were  measured.  To  measure 
the  lower  energy  neutrons,  the  experiment  was  moved  to  17.  59  meters 
and  the  electronics  gated  so  that  neutrons  from  327  keV  to  . 6 eV  were 
measured. 

Tables  14  and  15  present  the  spectra  of  source  neutrons  at  53.095 
and  17.  59  meters  respectively.  The  17.  59  meter  spectrum  also  includes 
the  high  energy  component  that  could  have  an  effect  on  the  results  by 


TABLE  14 


SPECTRUM  OF  NEUTRONS  AT  53.095  METERS 
(Supplied  by  IRT) 


TOTAL 

NEUTRON 


NEUTRON 

FLUX  IN  AND 

NEUTRON  FLUX 

FLUX  IN  BIN 

ABOVE  BIN 

(NEUTRONS/ 

(NEUTRONS/ 

(NEUTRONS/ 

LEADING 

AVERAGE 

MEV-SQ  CM- 

FRACTIONAL 

SQ  CM- 

SQ  CM- 

EDGE  ENERGY 

ENERGY 

SOURCE 

STD.  DEV. 

SOURCE 

SOUPCE 

(MEV) 

(MEV) 

MONITOR) 

OF  THE  FLUX 

MONITOR) 

MONITOR) 

1. 9872*01 

1.9163*01 

1.3930*05 

5.3127-02 

1.9764*05 

1. 9764*05 

1. 8454*01 

1.7818*01 

1. 5876+05 

5.0534-02 

2.0176*05 

3.9940*05 

1.7183*01 

1.6611*01 

2.0899*05 

4.6256-02 

2.3885+05 

6. 3825*05 

1. 6040*0 1 

1. 5524*01 

2.6479*05 

4.3404-02 

2.7311*05 

9.  1136*05 

1. 5008*01 

1. 4541*01 

3. 3305*05 

3.9794-02 

3.  1127*05 

1.2226+06 

1. 4074*01 

1.3649*01 

4.  1702*05 

3.7303-02 

3.5421*05 

1. 5768+06 

1. 3224*01 

1. 2837*01 

5.0042*05 

3.  4771-02 

3.8745*05 

1 . 9643+06 

1. 2450*01 

1.2096*01 

5.4763*05 

3.  2969-02 

3.8762*05 

2.3519+06 

1. 1742*01 

1. 1418*01 

6. 3690*05 

3.2105-02 

4.  1321*05 

2. 7651*06 

1. 1094*01 

1.0795*01 

6.7220*05 

3.  1714-02 

4.0076*05 

3.  1659*06 

1.0497*01 

1. 0223*01 

7.6227*05 

3.1252-02 

4.1860*05 

3. 5845+06 

9. 9482*00 

9.6947*00 

8.0478*05 

3.1220-02 

4.0797*05 

3. 9925*06 

9.4413*00 

9. 2066*00 

9.7173*05 

3.1041-02 

4.5574*05 

4.4482+06 

8. 9723*00 

8. 7549*00 

1.0145*06 

3.1922-02 

4.4104*05 

4.8892*06 

8. 5375*00 

8.3356*00 

1.  1442*06 

3. 2229-02 

4.6199*05 

5.  3512*06 

8. 1338*00 

7.9459*00 

1.  1848*06 

3. 3039-02 

4.4511*05 

5.7963*06 

7. 7581*00 

7. 5930*00 

1. 1725*06 

3.4250-02 

V , * (IQJ./JC 

6. 2069+06 

7.4079*00 

7.2445*00 

1.2865*06 

3.4806-02 

4.2052*05 

6.6274*06 

7.0811*00 

6. 9282*00 

1.3974*06 

3.5718-02 

4.2710*05 

7.0545*06 

6.7754*00 

6.6323*00 

1.3596*06 

3.6770-02 

3.8910*05 

7.4436*06 

6.4892*00 

6. 3550*00 

1. 5300*06 

3.  8770-02 

4.1064*05 

7.8542*06 

6. 2208*00 

6. 0948*00 

1.6221*06 

3.9680-02 

4.0879*05 

8. 2630*06 

5. 9688*00 

5.8504*00 

1.6110*06 

4.  1055-02 

3.8175*05 

8. 6448*06 

5. 7319*00 

5.6203*00 

1.6579*06 

4. 1821-02 

3.6986*05 

9.0146*06 

5. 5088*00 

5. 4036*00 

1.7023*06 

4.3996-02 

3. 5795*05 

9.3726*06 

5. 2986*00 

5.  1993*00 

1.8061*06 

4. 5086-02 

3.5839*05 

9.7310*06 

5.  1001*00 

5. 0063*00 

1. 8834*06 

4.5711-02 

3. 5305*05 

1 . 0084*07 

4.9126*00 

4.8240*00 

2.0608*06 

4.6697-02 

3.6536*05 

1.0449*07 

4. 7353*00 

4.6514*00 

2.3433*06 

4.7903-02 

3.9330+05 

1.0843+07 

4.5675*00 

4.4880*00 

2. 2979*06 

4.8602-02 

3.6551*05 

1.  1208*07 

4. 4094*00 

4.3330*00 

2.  3630*06 

4. 9726-02 

3.  5648*05 

1 . 1565*07 

4. 2576*00 

4. 1859*00 

2.4443*06 

4. 9779-02 

3.5013*05 

1.  1915*07 

4.  1143*00 

4. 0463*00 

2.4765*06 

4.9548-02 

3.  3709*05 

1. 2252*07 

3. 9782*00 

3. 9135*00 

2.5842*06 

5.0700-02 

3.  3456*05 

1.2586*07 

3. 8487*00 

3. 7871*00 

2.6925*06 

5.2077-02 

3. 3177*05 

1. 2918*07 

3. 7255*00 

3.6668*00 

2.7577*06 

5.2575-02 

3.2374*05 

1. 3242*07 

3. 6081 *00 

3. 5522*00 

2. 8068*06 

5.  2786-02 

3.  14  14*05 

1 . 3556*07 

3. 4962*00 

3.442e*00 

2.7090*06 

5.2491-02 

2. 9778*05 

1. 3854*07 

?. 3894*00 

3. 3385*00 

2.9541*06 

5. 5077-02 

3.01  18*05 

1.4155*07 

2875*00 

3.2388*00 

2.630  3*06 

5. 4363-02 

2.7576*05 

1 . 4431*07 

3.  1901*00 

3.  1435*00 

2.9004*06 

5.3601-02 

2.7012*05 

1. 4701*07 

3.0969*00 

3.0524*00 

2.9852*06 

5.  4906-02 

2.6603*05 

1.4967*07 

3. 0078*00 

2. 9651*00 

2.9931*06 

5.5699-02 

2.5537*05 

1. 5222*07 

2.9225*00 

2. 8816*00 

2.8698*06 

5.6905-02 

2.  3457*05 

1.5457*07 

TABLE  14.  (Cont'd.) 


TOTAL 

NEUTRON 


NEUTRON 

FLUX  IN  AND 

NEUTRON  FLUX 

FLUX  IN  BIN 

ABOVE  BIN 

(NEUTRONS/ 

(NEUTRONS/ 

(NEUTRONS/ 

LEADING 

AVERAGE 

MEV-SQ  CM- 

FRACTIONAL 

SQ  CM- 

SQ  CM- 

EDGE  ENERGY 

ENERGY 

SOURCE 

STD.  DEV. 

SOURCE 

SOURCE 

(MEV) 

(MEV) 

MONITOR) 

OF  THE  FLUX 

MONITOR) 

MONITOR) 

2. 8407400 

2.8016*00 

2.8767+06 

5.  809  1-02 

2.2538*05 

1.5682*07 

2. 7624400 

2. 7248*00 

3.0244+06 

5.6689-02 

2.2726*05 

1. 5910+07 

2. 6873400 

2.6512*00 

3.990  1+06 

5.  0786-02 

2.8776*05 

1 .6197*07 

2. 61 51400 

2. 5805*00 

4.6582+06 

4.6911-02 

3. 2257+05 

1. 6520+07 

2. 5459400 

2. 5126*00 

5. 1730*06 

4.6376-02 

3. 4409+05 

1 . 6864*07 

2. 4794400 

2.4474*00 

5.  1148  + 06 

4.4903-02 

3.2719+05 

1. 7191+07 

2. 41 54400 

2.3847*00 

5. 5766+06 

4. 3827-02 

3.4291*05 

1 .7534*07 

2. 3539400 

2. 3243*00 

5.6845*06 

4.3933-02 

3.  3649*05 

1.7871*07 

2. 2947400 

2. 2662*00 

6.4712+06 

4.2894-02 

3.6866*05 

1. 8239*07 

2. 2378400 

2. 2103*00 

6. 8677+06 

4. 2100-02 

3.7692*05 

1.8616*07 

2. 1829400 

2. 1564*00 

6.9771+06 

4.  1589-02 

3.6895*05 

1.8985+07 

2. 1300400 

2.0515*00 

6.7023+06 

4.0101-02 

1.0515*06 

2.0037+07 

*********** 

***  ABOVE 

DATA  TAKEN  WITH  NE-213  DETECTOR  ************** 

1. 9731400 

1.8866+00 

8.7712+06 

1. 5188-02 

1.5165*06 

2.  1553+07 

1. 8002400 

1.7246+00 

1 .0256+07 

1.4823-02 

1.5497*06 

2. 3103*07 

1.6491400 

1.5827*00 

1 . 1758+07 

1. 4554-02 

1.5615*06 

2.4664*07 

1. 51 63400 

1 . 4 57 6 + AO 

1 3074+07 

1. 4216-02 

1.5342+06 

2.6198+07 

1. 3990*00 

1.3468+00 

1 . 34c6  + 07 

1.4357-02 

1.4076+06 

2. 7606+07 

1. 2947*00 

1. 2482+00 

1.4171+07 

1. 4063-02 

1.3365*06 

2.8943+07 

1. 2017*00 

1 . 1600+00 

1.5251+07 

1 .3807-02 

1.2704*06 

3.0213*07 

1. 1184*00 

1 .0809+00 

1.5743*07 

1 . 3735-02 

1.  1807*06 

3.  1394  + 07 

1.0434*00 

1 . 0088  + 00 

1.6  176*07 

1. 3635-02 

1.  1181*06 

3. 2512*07 

9. 7428-01 

9. 4373-01 

1.7006*07 

1.3997-02 

1.0391*06 

3. 3551*07 

9. 13  18-01 

8.854 1-01 

1.9681*07 

1. 3864-02 

1.0929*06 

3. 4644*07 

8.5765-01 

8.3234-01 

2.0423*07 

1.4397-02 

1.0336*06 

3. 5677*07 

8. 0704:01 

7.8391-01 

2.0256*07 

1. 4907-02 

9.3704*05 

3.6614*07 

7.6078-01 

7. 3958-01 

1.9660*07 

1.5293-02 

8.3339*05 

3. 7448*07 

7. 1839-01 

6.9891-01 

2.  1250*07 

1.4817-02 

8.2769*05 

3.8275*07 

6. 7944-01 

6.6151-01 

2.0165*07 

1.5235-02 

7.2312*05 

3. 8999*07 

6.4358-01 

6.  2703-0  1 

1. 7720*07 

1.6265-02 

5.8653*05 

3. 9585*07 

6.  1048-01 

5. 9517-01 

1.4309*07 

1.8200-02 

4.3800*05 

4. 0023*07 

5.7987-01 

5.5252-01 

1. 6988*07 

1.  1889-02 

9.2907*05 

4. 0952+07 

5. 2518-01 

5.0152-01 

1.7795*07 

1.1572-02 

8.4188+05 

4.  1794*07 

4.7787-01 

4. 5727-01 

1. 8058*07 

1.  1258-02 

7.4381+05 

4. 2538*07 

4. 3668-01 

4.  1864-01 

1.9128*07 

1.1277-02 

6.90  14*05 

4. 3228*07 

4.0060-01 

3.8470-01 

2.0979*07 

1.  1044-02 

6.6692*05 

4.3895*07 

3. 6881-01 

3. 5473-01 

2. 1531*07 

1.  1 139-02 

6.063  1*05 

4. 4501*07 

3. 4065-*0  1 

3. 2813-01 

2.2053+07 

1.  1237-02 

5.5221*05 

4.5053*07 

3.  1561-01 

3.  044  1-01 

2. 1237+07 

1.1675-02 

4.7550*05 

4. 5529*07 

2. 9322-01 

2.8318-01 

2.0234+07 

1. 2255-02 

4.0630*05 

6.5935*07 

2. 7314-01 

2.6409-01 

2.3913*07 

1.  1652-02 

4.3259*05 

4 . 63  66  + 0’. 

2.5505-01 

2.4687-01 

2. 4 1 12*07 

1.2123-02 

3.9423*05 

4.6762+ 07 

2. 3870-01 

2.3129-01 

2.4621+07 

1.2691-02 

3.6488*05 

4.7127*07 

2. 2388-01 

2.  1713-01 

2. 2373*07 

1.  4229-02 

3.0181*05 

4.7429+07 

2.  1039-0  1 

2.0424-01 

2. 1205*07 

1.5776-02 

2.6082*05 

4.7690*07 
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TABLE  14.  (Cont'd.) 


NEUTRON  FLUX 
(NEUTRONS/ 


LEADING 

AVERAGE 

MEV- SQ  CM- 

EDGE  ENERGY 

ENERGY 

SOURCE 

(MEV) 

(MEV) 

MONITOR) 

1. 9809-01 

1.9246-01 

2. 2219*07 

1.0683-01 

1.8167-01 

2.4007*07 

1.7651-01 

1.7176-01 

2.4411*07 

1. 6702-01 

1.6264-01 

2.2955*07 

1. 5827-01 

1. 5423-01 

1.8849*07 

1.  5020-01 

1.4646-01 

2.0099*07 

1. 4273-01 

1. 3762-01 

2. 1 140*07 

1. 3252-01 

1.2794-01 

2.7561*07 

1.2337-01 

1. 1925-01 

2.7190*07 

1.  1514-01 

1. 1142-01 

2.8348*07 

1.0770-01 

1.0433-01 

2.790 1*07 

1. 0096-01 

9. 7900-02 

2.6380*07 

9.4840-02 

9. 2048-02 

2. 2712*07 

8. 9256-02 

8.6  7X13-02 

1.8576*07 

8.4151-02 

8.  1811-02 

2. 0088*07 

7. 9472-02 

7. 7322-02 

2.6703*07 

7. 5172-02 

7. 3192-02 

3. 2039*07 

7.  1213-02 

6. 9535-02 

3.5354*07 

6. 7558-02 

t>.  5867-02 

3.4131*07 

6. 4177-02 

6.2114-02 

3.6367*07 

NEUTRON 

TOTAL 
NEUTRON 
FLUX  IN  AND 

FLUX  IN  BIN 

ABOVE  BIN 

(NEUTRONS/ 

(NEUTRONS/ 

FRACTIONAL 

30  CM- 

SQ  CM- 

STD.  DEV. 

SOURCE 

SOURCE 

OF  THE  FLUX 

MONITOR) 

MONITOF) 

1.6541-02 

2.5019*05 

4.7940*07 

1.7031-02 

2.4775*05 

4. 8188*07 

1.7984-02 

2.3166*05 

4.8419*07 

1.9669-02 

2.0086*05 

4. 8620*07 

2. 2857-02 

1.5211*05 

4.8772*07 

2.3157-02 

1.5014*05 

4. 8922*07 

1. 9436-02 

2.  1584*05 

4.9138*07 

1.7746-02 

2.5218*05 

4. 9390*07 

1. 8670-02 

2.2377*05 

4.9614*07 

1.8900-02 

2.  1091  + 05 

4. 9825*07 

1.9651-02 

1.8805*05 

5.0013*07 

2.0850-02 

1. 6145*05 

5.0175*07 

2. 3310-02 

1.2682+05 

5.0301*07 

2.6671-02 

9.4830*04 

5.0396*07 

2. 6398-02 

9.3992*04 

5.0490+07 

2. 3529-02 

1 . 1482*05 

5.0605+07 

2. 2109-02 

1 . 2684*05 

5.0732+07 

2. 1631  C.7 

. ip. 2^05 

5.0861+07 

2. 2712-02 

1.  1540*05 

5.0976+07 

1.9728-02 

1. 5005*05 

5.  1126*07 
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TABLE  15.  SPECTRUM  OF  NEUTRONS  AT  17.  59  METERS 
(Supplied  by  IRT) 


TOTAL 

NEUTRON 

NF.U'i'Wl  FLUX 


FLUX 

NVj  i.I on 

NEUTRON 

IN  AND 

(NEUT/ 

FLUX 

FRACT- 

FLUX IN 

ABOVE 

UPPER 

LEADING 

MICRO- 

(NEUT/ 

IONAL 

BIN 

BIN 

EDGE 

EDGE 

AVERAGE 

AVERAGE 

SECOND- 

EV- 

STD. 

(NEUT/ 

(NEUT/ 

NEUTRON 

TIME 

NEUTRON 

TIME 

SQ  CM- 

SQ  CM- 

DEV. 

SQ  CM- 

SQ  CM- 

ENERGY 

(MICRO- 

ENERGY 

(MICRO- 

SOURCE 

SOURCE 

OF  THE 

SOUPCE 

SOUPCE 

(EV) 

SECONDS) 

(EV) 

SECONDS) 

MONITOR) MONITOR) 

FLUX 

MONITOR) MONITOR) 

1.9900*07 

2.8969-01 

1.8303*07 

3.0555-01 

3.85+07 

3.82-0 1 

8. 68-02 

1.  22*06 

1.22*06 

1.6707*07 

3.2142-01 

1.4685+07 

3.4123-01 

6.05+07 

5.94-01 

6.09-02 

2.40*06 

3.61*06 

1. 2663*07 

3.6103-01 

1.1451+07 

3.8084-01 

7.42*07 

1 . 21*00 

5.01-02 

2. 94+06 

6.55+06 

1.0240*07 

4.0065-01 

9.  1971*06 

4.2447-0 1 

8.52*07 

1. 94*00 

4. 39-02 

4.  05*06 

1.06*07 

3.  1543  *06 

4.4830-01 

7.4012*06 

4.7209-0 1 

8.  37*07 

2. 64*00 

4.73-02 

3.97*06 

1.46*07 

6. 64  82*06 

4.9587-01 

6.0069*06 

5.2365-01 

7.63*07 

3. 30*00 

5.04-02 

4.  24  + 06 

1.89*07 

5.3657*06 

5.5143-01 

4.8369*06 

5.8311-01 

6.  45*07 

3. 86*00 

5. 50-02 

4.09*06 

2.30+07 

4. 3081*06 

6. 147  9-01 

3.8809*06 

6.5052-0 1 

5.95*07 

4 . 97*00 

5.65-02 

4.25*06 

2.72+07 

3.4537*06 

6.8624-01 

3.1128*06 

7.2591-01 

4.78+07 

5. 53*00 

5.75-02 

3. 79*06 

3.09*07 

2.7720*06 

7.6557-01 

2.5019*06 

8.0918-01 

5.74*07 

9. 27*00 

4.73-02 

5.02*06 

3.59*07 

2.2319*06 

8.5278-01 

2.01 84*06 

9.0038-01 

7.63*07 

1.71*01 

3. 77-02 

7.  26*06 

4.32*07 

1 . 8050*06 

9.4798-01 

1.6201*06 

1 .0047  + 00 

8.  18*07 

2.51*01 

3. 37-02 

9. 27*06 

5.22*07 

k**»»***«* 

ABOVE  LATA  TAKEN  K 

ITH  NE-21 

2 DETECTOR  ***< 

>**«****! 

»**«**»* 

1.4353*06 

1.06  14*00 

1. 2899*06 

1. 1251+00 

7.  87*07 

3.45*01 

3.09-02 

1.00+07 

6.24*07 

1.1445*06 

1.1886*00 

1.0105+06 

1.2735*00 

6.88+07 

4. 35*01 

3.09-02 

1.  17*07 

7.40*07 

8. 7654*05 

1. 3583*00 

7.8464*05 

1.4430*00 

5.64*07 

5.22*01 

3. 09-02 

9. 59*06 

8. 37*07 

6.9274*05 

1.5279*00 

6. 2698*05 

1.6126*00 

3.62*07 

4. 66*0 1 

3.09-02 

6.  14*06 

8.98*07 

5. 61 22*05 

1.6975*00 

5.0236*05 

1 . 8034*00 

3.  11*07 

5.63*01 

3.09-02 

6. 60*06 

9.64*07 

4.4351*05 

1.9095*00 

4. 01 40*05 

2.0154*00 

2.79*07 

7. 04*0  1 

3.09-02 

5. 94*06 

1.02*08 

3.5930*05 

2.1215*00 

3.2288*05 

2. 2487*00 

1. 88*07 

6. 55*01 

3. 09-02 

4.76*06 

1.07*08 

2. 8647*05 

2.3758*00 

2.6010*05 

2.5031*00 

1.52*07 

7. 35*0  1 

3.09-02 

3. 86*06 

1. 11*08 

2.3373*05 

2.6303*00 

2.  1123*05 

2.7786*00 

1. 08*07 

7.  11  + 01 

3.09-02 

3. 20*06 

1.  14*08 

1 .8874*05 

2.9271*00 

1.70  15*05 

3.0967*00 

7. 63*06 

6. 94*01 

3.09-02 

2.59+06 

1.17*08 

1.5157*05 

3.2663*00 

1 . 3654*05 

3.4571*00 

5.36*06 

6. 80*01 

3. 09-02 

2. 05*06 

1.19*08 

1.2152*05 

3.6479*00 

1.0952*05 

3. 8599*00 

4. 28*06 

7. 55*0  1 

3.  18-02 

1.81*06 

1.21*08 

9.7527*04 

4.0719*00 

8.8019*04 

4. 3052*00 

2. 61*06 

6. 39*01 

3.62-02 

1.22*06 

1.22*08 

7.8511*04 

4. 5384*00 

7.0995*04 

4.7927*00 

2.34*06 

7.91*01 

3. 49-02 

1.  19*06 

1.23*08 

6.3479*04 

5.0472*00 

5.7536*04 

5. 3228*00 

2. 31*06 

1.07*02 

3. 24-02 

1. 27*06 

1.24*08 

5. 1594*04 

5.5984*00 

4.6885*04 

5.8952*00 

1.68*06 

1.06*02 

3. 55-02 

1.00*06 

1.25*08 

4.2176*04 

6.1920*00 

3.8217*04 

6.5312*00 

1.08*06 

9.23*0  1 

3.49-02 

7. 30+05 

1. 26*08 

3.4258*04 

6.8704+00 

3. 1004*04 

7.2520*00 

1. 10*06 

1.29*02 

3.65-02 

8. 42*05 

1.27*08 

2. 7750*04 

7.6336*00 

2.51  14*04 

8.0577*00 

1 . 10*06 

1.77*02 

3. 37-02 

9. 39*05 

1. 28*08 

2.2478*04 

8.4817*00 

2.0361*  J4 

8. 9481*00 

9.67*05 

2.  12*02 

3. 37-02 

8.98*05 

1. 29*06 

1. 8244*04 

9.4146*00 

1.6551*04 

9. 9234*00 

7. 42*05 

2. 23*02 

3.58-02 

7.55*05 

1.29*08 

1.4859*04 

1.0432*01 

1.3462*04 

1 .1005*01 

6.24*05 

2. 56*02 

3.62-02 

7.  16*05 

1.30*08 

1.2065*04 

1.1577*01 

1.0962*04 

1. 2191*01 

5. 31*05 

2.95*02 

3.68-02 

6.55*05 

1.31*08 

9.8595*03 

1.2807*01 

8.9602*03 

1.3485*01 

4.23*05 

3.  18*02 

3. 86-02 

5. 73*05 

1.31*08 

8.0609*03 

1.4164*01 

7. 3149*03 

1.4927*01 

3.67*05 

3. 76*02 

3. 83-02 

5. 58*05 

1.32*08 

6.5688*03 

1.5690*01 

5.9592+03 

1.6538*01 

2.92*05 

4.05*02 

3.99-02 

4.96*05 

1 .32*08 

5. 3497*03 

1.7386*01 

4.8564*03 

1 . e319*01 

2.58*05 

4. 89*02 

3. 99-02 

4. 81*05 

1.33*08 
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TABLE  15 


(Cont'd. ) 


UPPEP. 

EDGE 

LEADING 

EDGE 

AVEPAGE 

AVERAGE 

NEUTRON 

FLUX 

(NEUT/ 

MICPO- 

SECOND- 

NEUTRON 

FLUX 

(NEUT/ 

EV- 

FRACT- 

IONAL 

STD. 

NEUTRON 
FLUX  IN 
BIN 
(NEUT/ 

TOTAL 
NEUTRON 
FLUX 
IN  AND 
ABOVE 
BIN 
(NEUT/ 

NEUTRON 

TIME 

NEUTRON 

TIME 

so  CM- 

SQ  CM- 

DEV. 

SQ  CM- 

SQ  CM- 

ENERGY 

(MICRO- 

ENERGY 

(MICRO- 

SOURCE 

SOURCE 

OF  THE 

SOURCE 

SOURCE 

(EV) 

SECONDS) 

(EV) 

SECONDS) 

MONITOR) MONITOR) 

FLUX 

MONITOR) MONITOR) 

4.3630*03 

1.9251*01 

3.9658*03 

2.0269*01 

2.20*05 

5.63*02 

4.05-02 

4. 48*05 

1.33+08 

3. 5686*03 

2.1287*0 1 

3.2389*03 

2.  2432  + 0 1 

1 . 84*05 

6.  39*02 

4.  1 1-02 

4.20*05 

1.34*08 

2.9091*03 

2. 3576*01 

2. 6435*03 

2.4827*01 

1 . 54*05 

7. 26*02 

4. 23-02 

3. 86*05 

1.34*08 

2.3778*03 

2.6076*01 

2. 1500*03 

2.7541*01 

1 . 32*05 

8. 47*02 

4.29-02 

3.  86*05 

1. 35*08 

1.9223*03 

2. 9003*01 

1.72  16*03 

3.0806*01 

1 . 09*05 

9. 73*02 

4.20-02 

3.92*05 

1.35*08 

1. 5208*03 

3.2608*01 

1.3770*03 

3.4410*01 

9.46*04 

1.  18*03 

4.42-02 

3.  40+05 

1 .35*08 

1.2332*03 

3.6212*01 

1 . 1086*03 

3.8375*01 

8.03*04 

1.  39*03 

4.33-02 

3. 46*05 

1.36+08 

9.8405*02 

4.0538*01 

8.9376*02 

4. 2700*01 

6.64*04 

1 . 59*03 

4.70-02 

2.87*05 

1.36*08 

8.0346*02 

4.4863*01 

7.2633*02 

4.7385*01 

5.74*04 

1.89*03 

4. 60-02 

2.  91*05 

1. 36*08 

6.4921*02 

4. 9908*01 

5.  8545*02 

5.2792*01 

5.  10*04 

2.31*03 

4. 51-02 

2.95*05 

1. 36*08 

5.2169*02 

5.5674*01 

4.7008*02 

5. 89  18*0  1 

U.  19^04 

2.63*03 

4.70-02 

2.71+05 

1 .37*  08 

4. 1848*02 

6.2162*01 

3.7726*02 

6.5766*01 

3.70*04 

3. 25*03 

3.18-02 

2. 68*05 

1. 37+08 

3.3603*02 

6.9370*01 

3.0333*02 

7.3335+01 

3. 35*04 

4.07*03 

3.09-02 

2.66*05 

1.37*08 

2.7063*02 

7. 7299*01 

2.4476*02 

8.  1624*0  1 

2. 85*04 

4.76*03 

3.09-02 

2. 46*05 

1.37*03 

2. 1890*02 

8. 5949+01 

1.9844*02 

9.0634*0 1 

2. 49*04 

5. 68*03 

3.09-02 

2.  33*05 

1. 38*08 

1.7748*02 

9.5319*01 

1 . 6 07  ^ *0  2 

1. 0072+02 

•» . *2-  * 0 4 

6.14*03 

3. C9-C2 

2.10-05 

1 . 3 3+  G3 

1.4356*02 

1.06  13*02 

1.3018*02 

1 . 1189*02 

"i  . o 0 ♦ 0 4 

6.89*03 

i.  24  -02 

1 . 84*05 

1.39*08 

1.1680*02 

1. 1767*02 

1.0577+02 

1 . 2416*02 

1.47+04 

8. 62*03 

3.  12-02 

1. 90*05 

1. 39*08 

9.4748*01 

1. 3064*02 

8.5800*01 

1. 3785+02 

1. 23*04 

9.90*03 

3.  12-02 

1 . 77*05 

1.39*08 

7.6851*01 

1.4506*02 

6.9651*01 

1.5298*02 

1 . 04*04 

1 . 15*04 

3.  15-02 

1. 66*05 

1 .39*08 

6.2450*01 

1.6091*02 

5.6682*01 

1 .6956  + 02 

9.86*03 

1. 48*04 

3. 09-02 

1.71*05 

1. 39*08 

5.0914*01 

1.7821*02 

4.6148*01 

1.  8795*02 

8. 23*03 

1 . 67*04 

3.09-02 

1 . 60*05 

1.39*08 

4. 1382*01 

1.9767*02 

3.76  14*01 

2.0812+02 

7.54*03 

2. 08*04 

3.09-02 

1. 57*05 

1.40*06 

3.3846*01 

2. 1858*02 

3.0687*01 

2. 3047*02 

5. 45*03 

2.05*04 

3.43-02 

1.  30*05 

1.40*08 

2.7528*01 

2.4236*02 

2.5538*01 

2.5221*02 

5. 00*03 

2.48*04 

3.37-02 

9.90*04 

1.40*  08 

1 ABOVE  DATA  TAKEN  WITH  LI-6  GLASS  SCINTILLATOR  DETECTOR  ***< 

k *****  ** 

2.3547*01 

2. 6206*02 

2.1290*01 

2.767-7*02 

3.24*03 

2.  12*04 

4.51-02 

9. 54*04 

1.40*08 

1.9034*01 

2.9147*02 

1.7240*01 

3.0752*02 

3.50*03 

3.  14*04 

3. 92-02 

1.  12*05 

1 . 40*03 

1.5446*01 

3.2356*02 

1.4020*01 

3. 4094*02 

3.01*03 

3. 68*04 

3.86-02 

1.05*05 

1.40*08 

1.2594*01 

3.5833*02 

1.1390*01 

3.7838*02 

3.24*03 

5.43*04 

3.  34-02 

1.  30*05 

1 . 40*08 

1.0186*01 

3.9844*02 

9.2462*00 

4.1983*02 

2. 76*03 

6. 29*04 

3. 37-02 

1.  18*05 

1.40*  08 

8.3064*00 

4.4122*02 

7. 5295*00 

4.6529*02 

2.  12*03 

6. 55*04 

3.46-02 

1.02*05 

1.40*08 

6.7527*00 

4.8936*02 

6.1202*00 

5. 1609*02 

7. 28*02 

3.09+04 

5.59-02 

3. 89*04 

1.41*08 

5.4877*00 

5.4283*02 

4.9773*00 

5.7225*02 

1. 51*03 

8.72*04 

3.43-02 

8.38*04 

1.41*08 

4.4670*00 

6.0167*02 

4 .0572*00 

6. 3375*02 

1.42*03 

1.11*05 

3.21-02 

9.08*04 

1.41*08 

3.6474*00 

6.6584*02 

3.3081*00 

7. 01 94*02 

1 . 24*03 

1. 32*05 

3.09-02 

8. 93*04 

1. 4 1*08 

2.9687*00 

7.3804*02 

2.6919*00 

7.7815*02 

1.  11*03 

1.62*05 

3.00-02 

8. 98*04 

1.4  1+  00 

2.4151*00 

8.  1 827*02 

2.1914*00 

8.6239*02 

9. 67*02 

1 . 90*05 

3. 03-02 

8.52*04 

1.41*08 

1.9678*00 

9.0650*02 

1.7880*00 

9.5465*02 

8.37*02 

2. 23*05 

3.00-02 

8.06*04 

1.41*08 

1.6081*00 

1.0028*03 

1.4606*00 

1.0563*03 

7.  18*02 

2.61*05 

2.97-02 

7.67*04 

1.41*08 

1. 31  31*00 

1.  1097*03 

1 . 1 933*00 

1 .1686*03 

6. 19*02 

3.04*05 

2.97-02 

7. 30*04 

1.41*08 

1.0734*00 

1.2274*03 

9.7485-01 

1.2929*03 

5.31*02 

3. 53*05 

3.00-02 

6. 94*04 

1. 41*08 

8.7630-01 

1. 3585*03 

7.96  14-0  1 

1. 4306*03 

4. 36*02 

3.92*05 

3. 06-02 

6. 29*04 

1.41*08 
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TABLE  15 


(Cont'd. ) 


UPPER 

EDGE 

LEADING 

EDGE 

AVERAGE 

AVERAGE 

NEUTRON 

FLUX 

(NEUT/ 

MICRO- 

SECOND- 

NEUTRON 

FLUX 

(NEUT/ 

EV- 

fract- 

ional 

STD. 

NEUTRON 
FLUX  IN 
BIN 
(NEUT/ 

TOTAL 
NEUTRON 
FLUX 
IN  AND 
ABOVE 
BIN 
(NEUT/ 

NEUTRON 

TIME 

NEUTRON 

TIME 

SQ  CM- 

SQ  CM- 

DEV. 

SO  CM- 

SO  CM- 

ENERGY 

(MICRO- 

ENERGY 

(MICRO- 

SOURCE 

SOURCE 

OF  THE 

SOUPCE 

SOUPCE 

(EV) 

SECONDS) 

(EV) 

SECONDS) 

MONITOR)  MONITOR) 

FLUX 

MONITOR) MONITOR) 

.1599-01 

1.5028*03 

6.5026-01 

1. 58  30  + 03 

3. 42*02 

4.  15*05 

3.28-02 

5. 48*04 

1.41+08 

i.  8452-0  1 

1. 6633*03 

5.3114-01 

1.7515*03 

2.55+02 

4.22  + 0 5 

3.62-02 

4.  50*04 

1.41*08 

.7776-01 

1.8398*03 

4.  3413-01 

1.9373*03 

1.55+02 

3.45+05 

4.60-02 

3.02*04 

1.4  1*08 

ABOVE  DATA  TAKEN  WITH  BF3  PROPORTIONAL  COUNTER  *************** 


thermalizing  in  the  experiment  and  residing  until  later  times.  For  this 
reason  our  calculations  at  17.53  meters  included  the  entire  spectrum 
and  results  were  also  calculated  at  times  earlier  than  the  experimental 
times.  This  source  data  was  supplied  by  IRT  on  computer  cards  which 
were  used  as  input  to  our  source  routines  (described  in  the  next  section^. 
The  beam  profile  measured  by  IRT  was  also  included  in  our  calculations. 

The  IRT  silicon  dosimeter  is  a silicon  surface  barrier  detector 
surrounded  by  silicon  and  encased  in  an  aluminum  can.  The  silicon  surface 
barrier  detector  measures  ionizing  energy  deposition  above  200  keV,  so 
that  the  entire  silicon  dosimeter  has  a response  much  like  a free  field 
ionizing  silicon  response  which  has  been  attenuated  through  the  surround- 
ing silicon  and  aluminum  can. 

The  response  of  the  dosimeter  to  neutrons  and  to  gamma  rays  was 
supplied  to  us  by  IRT.  The  response  to  neutrons  was  measured,  and  the 
response  to  gamma  rays  was  calculated  with  the  SANDYL  code  and  verified 
by  measurement.  Our  approach  was  to  use  the  IRT  response  at  the  out- 
side of  the  dosimeter  and  also  to  include  the  dosimeter  in  our  calculations. 
This  required  a special  treatment  of  the  dosimeter.  For  all  adjoint  cal- 
culations, the  adjoint  particles  were  initialized  isotropically  at  the  center 
of  the  dosimeter  sampling  the  IRT  response  to  determine  the  energy  and 
particle  type.  For  the  particle's  first  flight,  the  dosimeter  was  assumed 
void  so  that  the  dosimeter  response  would  be  maintained  to  the  surface  of 
the  dosimeter.  For  the  remainder  of  the  adjoint  random  walk,  the  dosi- 
meter material  was  included. 

For  the  forward  calculations,  which  used  a next  flight  estimator, 
the  detector  was  assumed  void  for  ail  estimations  and  the  material  was 
included  for  all  transport.  Also  no  estimation  was  performed  for  scatter- 
ings within  the  dosimeter.  These  treatments  assured  that  the  dosimeter 
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response  would  apply  to  the  surface  of  the  dosimeter  and  that  the  dosimeter 
be  included  in  the  transport  calculation. 

Our  dosimeter  response  was  presented  in  the  previous  section,  see 
Table  5,  in  our  36  energy  group. 
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4.  STEEL-STYROFOAM  CALCULATIONS 


Before  attempting  the  more  complex  Sprint  missile  experiment, 

IRT  fabricated  and  measured  the  steel-styrofoam  cylindrical  sample  which 
was  designed  to  have  characteristics  similar  to  the  real  Sprint  experiment. 
Using  the  same  procedures  later  used  for  the  Sprint  experiment,  the  cy- 
linder response  was  measured  at  both  the  17  meter  and  the  53  meter 
stations  of  the  flight  path.  The  17  meter  measurement  observed  the  re- 
sponse to  the  lower  energy  neutrons  from  . 25  MeV  to  . 4 eV,  and  the  53 
meter  experiment  measured  the  higher  count  rates  of  the  high  energy 
neutrons  from  100  keV  to  20  MeV.  In  this  section  details  and  results  of 
three  time  dependent  Monte  Carlo  calculations  of  these  two  experiments 
are  presented.  The  17  meter  experiment  was  calculated  with  both  forward 
and  adjoint  techniques  while  the  50  meter  experiment  was  calculated  only 
with  the  adjoint  method. 

4.  1 FORWARD  MONTE  CARLO  TECHNIQUES 

The  forward  Monte  Carlo  calculations  of  the  cylinder  experiment 
was  performed  with  the  standard  version  of  MORSE  which  was  modified 
slightly  to  treat  the  time  dependent  nature  of  the  source.  The  scoring 
technique  used  was  the  next  flight  estimator  (subroutine  RELCOL).  A 
track  length  per  volume  estimator  was  also  included  as  a check.  The 
LINAC  source  as  measured  by  IRT  was  previously  presented  in  Tables 
14  and  15.  Description  of  the  Monte  Carlo  techniques  for  both  the  forward 
and  adjoint  calculations  follow.  These  same  techniques  were  also  employed 
for  the  missile  calculations  which  are  described  in  the  next  section. 
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Our  actual  sampling  scheme  was  to  select  a continuous  energy  off 
a uniform  lethargy  distribution  and  adjust  the  source  weight  to  agree  with 
the  IRT  source.  The  energy  was  selected  from 


E 


= E . 

min 


where  £ is  a uniform  random  number  between  0 and  1 and  E . and 
s min 

E are  the  lowest  and  highest  energies  respectively.  The  source 
nicix 

weight  was  then  calculated  by 


SW  = A • E • f(E)  in 


max 

min 


which  assures  a fair  game.  Here  f(E)  is  the  Intelcom  source  in  units 
2 

neutrons/(eV-cm  -source  monitor)  and  A is  the  area  of  the  beam.  Rather 
than  setting  the  source  particle’s  weight  to  SW,  the  weight  was  set  to  1.0 
and  SW  was  carried  with  the  particle  and  all  its  secondaries  to  be  in- 
cluded later  as  a factor  in  all  scores.  In  this  manner  transport  importance 
parameters  were  independent  of  a particle's  initial  energy. 

The  actual  steps  taken  in  the  forward  MORSE  calculation  are  listed 

below: 

1.  A continuous  energy  E was  sampled  from  a uniform  lethargy 
distribution  and  a source  weight  was  calculated  as  a function 
of  E and  the  Intelcom  supplied  source. 

2.  The  time  of  arrival  t at  the  front  of  the  cylinder  was  cal- 
culated from  E and  the  flight  path  length. 

3.  The  energy  group  IG  corresponding  to  E was  determined. 
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4.  An  uncollided  contribution  was  scored  in  the  time  bin  cor- 
responding to  the  time  of  arrival  at  the  detector  center. 
The  score  was 


^ F(IG)  R(IG) 

where  the  F(IG)  are  pre -calculated  attenuation  factors  for 
each  energy  group  through  the  centerline  of  the  sample  to  the 
detector  point  and  R(IG)  is  the  detector  response  (Table  5). 

5.  A source  position  was  sampled  on  the  front  surface  of  the  cy- 
linder from  the  beam  shape. 

6.  From  here  a standard  MORSE  calculation  with  WATE  = 1.0 
and  using  the  next  flight  estimator  was  initiated  with  the  in- 
itial values  calculated  in  steps  1-5  above. 

7.  For  each  score  the  source  weight  factor  was  also  included. 

4.  2 TIME  DEPENDENT  ADJOINT  CALCULATION 

In  Section  2 the  adjoint  next  flight  estimator  to  a beam  was  dis- 
cussed for  the  calculation  of  the  steady  state  energy  response.  The  ex- 
tension of  this  technique  to  calculate  the  time  dependent  deposition  should 
be  straightforward,  especially  since  the  time  a particle  takes  to  travel 
from  its  source  to  the  detector  will  be  the  same  as  an  adjoint  history 
traveling  the  reverse  trajectory  from  the  detector  to  the  source.  In  a 
time  dependent  adjoint  calculation  the  next  flight  estimator  knows  a particle's 
weight  (WATE),  an  energy  group  (IG),  and  the  time  from  the  detector  (AGE) 
when  the  particle  is  entering  the  beam  at  the  front  of  the  cylinder.  It  only 
remains  to  determine  the  correct  time  bin  (or  time  bins)  and  source  in- 
tensity for  the  score.  A misconception  at  this  point  of  how  to  fold  the 
adjoint  fluence  with  the  time  and  energy  dependent  source  caused  consider- 
able consternation  and  as  such  is  worth  explaining  in  detail. 
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The  weight  of  the  adjoint  particle  upon  entering  the  source  beam 
represents  a total  dose  in  the  detector  per  incident  particle  per  unit  area 
from  a source  neturon  with  energy  group  IG  [the  units  of  weight  are  rad  (Si)/ 
(n/cm  )].  The  time  AGE  represents  a random  estimate  of  time  elapsed 
between  entry  of  the  neutron  into  the  sample  and  deposition  of  the  dose. 

The  three  plots  of  Fig.  5 schematically  depict  the  correct  procedure  for 
scoring  the  time  dependence. 

Plots  a and  b define  the  source;  a)  giving  the  energy -time  corre- 
lation which  in  this  case  is  defined  by  the  flight  path,  and  b)  depicting  the 
time  dependence  of  the  source.  In  our  case  of  energy  and  time  correlation 
these  two  curves  exactly  define  the  source.  More  complicated  energy  and 
time  dependent  sources  would  be  represented  by  a surface  in  time -energy - 
source  intensity  space.  The  adjoint  radiation  transport  results  are  con- 
tained in  the  three  quantities,  AGE  (transport  time),  WATE  (adjoint  particle 
weight),  and  IG  (the  energy  group  number).  The  third  plot  c)  then  defines 
the  correct  transformation  required  to  fold  the  transport  results  with  the 
time  dependent  source.  The  energy  group  IG  on  plot  a)  defines  a time 
domain  on  plot  b)  which  encloses  a portion  of  the  time  dependent  source 
rate.  It  is  this  portion  of  the  time  dependent  source  rate  that  should  be 
multiplied  by  WATE  and  scored  as  a dose  rate  at  the  times  AGE  later  than 
the  source  times.  Notice  that  the  shape  of  the  neutron  source  rate  should 
be  preserved  when  scoring  it  into  the  dose  rate  results. 

For  actual  computer  implementation  the  process  is  complicated 
by  the  existance  of  three  sets  of  time  bins  associated  with  each  of  the  three 
curves.  The  energy  group  boundaries  define  broad  time  boundaries  TIG 
corresponding  to  the  times  of  arrival  of  the  energy  group  edges.  The 
source  rate  actually  exists  as  a fine  histogram  with  times  TS0R  and  av- 
erage count  rates  RS0R.  The  final  count  rate  then  is  to  be  scored  into 
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Figure  5 


Graphical  Depiction  of  Adjoint  Time  Dependent  Scoring. 

(a)  Source  Energy-Time  Correlation,  (b)  Time  Dependence 
of  the  Source,  and  (c)  Procedure  for  Scoring  a Time 
Dependent  dose  rate 


time  bins  TBIN  decided  by  the  code  user.  Figure  6 depicts  these  three 
time  bins  and  demonstrates  our  algorithm  for  scoring  the  dose  rate. 
Before  any  transput  is  performed  the  source  rate  histogram  (TS0R, 
RS0R)  is  represented  by  TSAV(IG,  I)  and  S0R(IG,  I)  where  the  TSAV's 
are  the  centers  of  the  smallest  time  bins  bounded  either  by  the  TS0R’s 
or  the  TIG's.  The  S0R’s  then  are  the  total  number  of  neutrons  repre- 
sented at  times  TSAV  which  are  the  product  of  the  corresponding  RS0R 
and  time  product  of  the  corresponding  RS0R  and  time  bin  widths.  Let 
NIS(IG)  be  the  number  of  times  TSAV  that  lie  between  TIG(IG)  and  TIG 
(IG-1).  For  each  adjoint  estimate  to  the  beam,  NIS(IG)  scores  are  made 
at  times  AGE  + TSAV  (IG,  I).  Each  scoring  value  is  WATE*S0R(IG,  I). 
The  scoring  routine  determines  which  time  bin  TBIN  contains  AGE  +TSAV 
(IG,  I)  and  sums  the  scoring  value  in  the  appropriate  bin.  At  the  conclu- 
sion of  the  Monte  Carlo  calculation,  each  value  in  the  TBIN’s  array  is 
divided  by  the  time  bin  widths  TBIN (I)  -TBIN(I  + 1)  to  change  the  units 
back  to  a dose  rate  (rads/sec).  This  procedure  performs  the  transfor- 
mation indicated  in  the  three  curves  of  Fig.  5 within  the  resolution  of 
the  user  defined  time  bins  TBIN. 

When  first  attempting  the  time  dependent  adjoint  calculation  we 
unfortunately  discovered  that  other  more  intuitive  scoring  algorithms 
can  in  fact  given  the  wrong  answer.  It  is  essential  that  the  shape  of  the 
time  dependent  source  rate  possible  with  group  IG  be  scored  in  the  anal- 
ysis time  bins  with  the  same  shape  but  with  time  delayed  by  AGE.  We 
maintained  the  shape  by  scoring  several  times  for  each  estimate.  A 
single  time  could  be  sampled  and  the  entire  source  enclosed  by  TIG(IG) 
and  TIG(IG  + 1)  scored  at  this  time,  but  of  course  this  would  increase  the 
variance. 


71 


RESULT 

RADS/SEC 


TIG(IG)  TIG(IG-I) 


Figure  6.  Three  Sets  of  Time  Bins  Used  for  Adjoint 
Time  Dependent  Scoring 
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This  adjoint  technique  was  used  to  calculate  the  steel-styrofoam 
and  Sprint  missile  time  dependent  response  for  both  the  17  and  53  meter 
experiments. 

4.  3 COMPARISONS  OF  THE  STEEL-STYROFOAM  RESULTS  - 

FORWARD,  ADJOINT,  AND  EXPERIMENT 

The  steel -styrofoam  17  meter  experiment  was  calculated  with 
both  forward  and  adjoint  Monte  Carlo  while  the  53  meter  experiment  was 
calculated  with  the  adjoint  technique  only.  Figure  7 compares  the  adjoint 
steel-styrofoam  calculation  at  53.  59  meters  to  the  IRT  experiment.  The 
results  are  plotted  as  dose  with  the  units  rad(Si)/(psec  • source  monitor) 
versus  time  ( p sec  ).  The  verticle  bars  represent  a standard  deviation. 
The  53  meter  adjoint  calculation  tracked  40,000  histories  in  5.  77  minutes 
on  a CDC  7600  computer.  Splitting  and  Russian  roulette  parameters  were 
set  so  that  10%  of  the  particles  split  and  50%  were  killed  by  Russian  rou- 
lette. 

The  agreement  with  the  experiment  appears  to  be  excellent  up  to 
4 microseconds  where  the  calculation  drops  40%  below  the  experiment. 
Times  after  4 microseconds  corresponds  to  incident  neutron  energies  from 
1 to  . 1 MeV.  These  neutrons  contribute  small  amounts  of  neutron  ioniza- 
tion, almost  no  inelastic  gamma  rays,  and  must  loose  considerable  energy 
before  they  generate  capture  gamma  rays.  Therefore,  this  component 
of  the  dose  is  both  small  and  difficult  to  calculate.  The  excellent  agree- 
ment above  1 MeV  demonstrates  that  the  fast  neutron  effects  in  iron, 
polystyrene  and  silicon  are  well  understood  and  can  be  accurately  calcu- 
lated. 

The  17  meter  steel -styrofoam  experiment  was  calculated  with  both 
the  forward  and  adjoint  methods.  As  such,  it  represents  not  only  a good 
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Figure  7. 


Comparison  of  the  53  Meter  Steel -Styrofoam  Adjoint 
Calculation  With  the  IRT  Linac  Experiment 


compairson  between  experiment  and  calculation  but  also  an  excellent  com- 
parison of  two  calculational  techniques.  Both  calculations  used  very 
little  Russian  roulette  or  splitting  as  the  primary  history  termination 
was  escape.  For  the  forward  calculation  48,000  histories  were  tracked 
in  8.  73  minutes  on  a CDC  7600  computer  with  an  average  of  eleven  scat- 
terings per  history.  The  adjoint  calculation  tracked  80,000  histories  in 
12. 15  minutes  on  the  same  computer  with  splitting  and  Russian  roulette 
parameter  set  so  that  10%  of  the  particles  split  and  50%  were  killed  by 
Russian  roulette. 

Figure  8 compares  the  forward  and  adjoint  17  meter  calculations 
with  the  experiment.  Usually  comparisons  of  forward  and  adjoint  Monte 
Carlo  calculations  are  not  time  dependent.  As  such,  the  only  possible 
comparison  is  between  the  two  total  doses,  two  single  numbers  that  may 
be  dominated  by  special  aspects  of  the  transport  and  therefore  not  repre- 
sentative of  the  entire  calculations.  The  time  dependent  results  presented 
here  show  the  entire  character  of  the  two  calculations.  For  instance  dose 
before  1.3  psec  is  from  fast  neutrons  and  inelastic  gamma  rays  while 
the  dose  after  500  psec  is  due  to  thermal  capture.  Between  these  times 
epithermal  capture  and  thermal  capture  are  both  contributing  dose. 

Returning  to  the  forward  adjoint  comparison,  first  notice  that  the 
forward  calculation  is  generally  lower  than  the  adjoint  results,  at  spots 
as  much  as  50%.  This  probably  represents  the  nature  of  the  next  flight 
estimator  in  the  forward  calculation  which  tends  to  undersample  collisions 
close  to  the  dose  point.  In  this  steel-styrofoam  calculation,  most  of  the 
collision  density  occurs  in  the  steel  shell,  and  except  for  low  energy  neu- 
trons, there  is  little  scattering  in  the  styrofoam  or  the  dosimeter.  There- 
fore the  next  flight  estimator  generates  resonable  statistics  with  no  infinite 


75 


DOSE  RATE  <RAD(Si)/MICROSECOND  - SOURCE  MONITOR) 


Figure  8.  Comparison  of  the  17  Meter  Steel-Styrofoam  Forward 
and  Adjoint  Calculations  with  the  Linac  Experiment 


variance  signs,  but  may  still  be  underestimating  the  contribution  of  col- 
lisions near  the  dosimeter. 

It  is  also  of  interest  to  notice  the  relative  efficiencies  of  the  for- 
ward and  adjoint  calculations.  In  general  it  is  expected  that  a forward 
calculation  will  be  more  efficient  when  the  source  is  small  and  the  detector 
is  large,  and  an  adjoint  calculation  will  be  more  efficient  when  the  detector 
is  small  and  the  source  is  large.  Thus,  for  the  steel-styrofoam  calcula- 
tion, the  adjoint  technique  might  be  expected  to  be  more  efficient.  How- 
ever, Fig.  7 shows  that  at  times  around  250  psec  the  adjoint  fractional 
standard  deviation  is  0.  270  for  12.  15  minutes  of  computing  while  the  for- 
ward fsd  is  0.  128  for  8.  73  minutes  of  computing.  Defining  a relative 
efiiciency  as  tj(fsdj)  /t2(fsd2)  , the  adjoint  calculation  would  require 
6. 19  more  computing  time  to  calculate  the  same  variance  as  the  forward 
calculation.  The  explanation  of  this  startling  result  lies  in  the  estimates 
that  are  made  during  thermal  scattering.  In  the  forward  random  walk, 
estimates  of  capture  gamma  ray  dose  contributions  are  made  from  many 
of  the  thermal  neutron  collision  sites.  Since  thermal  collision  each 
occurs  at  a different  time,  estimates  are  made  to  several  of  the  time 
bins  from  a single  random  walk.  However,  in  the  adjoint  random  walk 
meaningful  estimates  cannot  be  made  from  the  thermal  collisions.  The 
upscatter  probabilities  are  such  that  the  next  flight  estimation  can  only 
be  made  to  the  energy  group  immediately  above  the  thermal  group.  Since 
there  is  no  source  in  either  this  group  or  the  thermal  group,  next  flight 
estimation  gives  no  sccre.  The  only  information  gained  from  the  thermal 
random  walk  is  the  location  and  time  of  the  eventual  upscatter  or  escape. 
Once  a single  upscatter  from  the  thermal  group  has  occurred,  the  time 
wiU  not  change  significantly  for  the  rest  of  the  history.  Thus  a single 
thermal  random  walk  (perhaps  50  collisions)  will  only  generate  a few 
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scores  in  a few  adjacent  time  bins  (or  no  score  at  ail  if  leakage  occurs 
before  upscattering).  Clearly  the  general  feeling  regarding  detector  sizes 
and  forward  versus  adjoint  calculation  efficiencies  is  not  always  true  for 
time  dependent  calculations.  In  this  case  the  ability  of  the  forward  next 
flight  estimator  to  generate  more  estimates  over  longer  times  makes  the 
forward  calculation  more  efficient  than  the  adjoint. 

Returning  to  Fig.  8,  notice  that  both  the  forward  and  the  adjoint 
calculations  are  lower  than  the  experiment.  In  general,  the  calculations 
are  30%  lower  than  the  experiment  but  at  a few  times  they  are  as  much  as 
70%  lower  than  the  experiment.  Also  the  spectral  shapes  agree  approxi- 
mately within  the  statistics  of  the  calculations. 

Possible  sources  of  this  difference  between  experiments  and  cal- 
culations could  result  from  any  of  the  following  effects: 

1.  Experimental  normalization 

2.  Too  much  absorption  or  leakage  in  the  calculations 

3.  Incorrect  detector  response  to  low  energy  neutrons 

4.  Incorrect  capture  gamma  ray  yields  or  energy  spectra 

5.  Other  calculational  difficulties 

Another  forward  calculation  with  slightly  less  absorption  and  a slightly 
harder  capture  gamma  ray  spectrum  showed  very  little  difference  from 
the  forward  calculation  of  Fig.  8,  therefore,  eliminating  reasons  2 and 
4.  Reasons  1 and  5 would  be  difficult  to  verify,  but  both  IRT  and  SAI 
used  considerable  care  in  our  procedures  to  eliminate  such  errors.  The 
normalization  of  the  experimental  results  depends  on  an  accurate  determin- 
ation of  the  dosimeter's  active  volume  which  did  prove  to  be  a difficult  ex- 
perimental task.  The  good  experiment -calculation  agreement  at  high 


energies  indicates  that  the  experimental  normalization  is  correct.  The 
difference  was  found  to  be  reason  3),  incorrect  detector  response  to  low 
energy  neutrons.  Unfortunately  this  fact  was  not  found  until  after  the 
Sprint  missile  calculations  were  completed  and  the  dosimeter  treatment 
already  described  was  used  in  all  the  missile  calculations. 

4.  4 IMPROVED  DOSIMETER  TREATMENT 

The  17  meter  calculations  of  the  steel-styrofoam  sample  (Fig.  8) 
and  the  Sprint  missile  indicated  that  the  dosimeter  treatment  described 
in  Section  3.  4 and  tabulated  in  Table  5 was  generating  low  results  at  late 
times  by  ignoring  capture  in  the  detector.  Originally  it  was  planned  to 
include  all  effects  of  the  detector  by  including  it  in  the  calculations  and 
using  a free  field  silicon  dose  response  at  the  detector’s  center.  This 
free  field  response  was  to  be  corrected  for  low  energy  gamma  ray  attenua- 
tion in  the  detector,  thus  allowing  the  detector  to  be  treated  as  a point. 

This  technique  proved  inaccurate  for  neutrons  from  0.5  to  2 MeV  where 
the  95  keV  detector  bias  prevented  accurate  determination  of  a neutron 
ionizing  response.  Because  of  this  we  elected  to  use  the  measured  dos- 
imeter response  at  all  times.  In  measuring  the  detector  response  to  low 
energy  neutrons  IRT  found  no  response  to  neutrons  below  300  keV  with 
the  single  exception  of  a gold  line  between  4.44  and  5.36  eV.  For  this 
reason  we  assumed  that  the  17  meter  results  would  be  dominated  by  cap- 
ture in  the  steel  and  the  styrofoam. 

To  test  the  effect  of  capture  in  the  detector,  a simple  modification 
to  a forward  17  meter  calculation  was  made.  As  discussed  in  Section 
4.  1,  scatterings  inside  the  detector  were  ignored  because  their  effect 
was  already  included  in  the  next  flight  estimator  contributions  which  as- 
sumed the  dosimeter  response.  The  modification  to  include  capture  in 
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the  detector  was  to  perform  the  next  flight  estimation  of  a capture  gamma 
ray  to  the  detector  center  for  each  neutron  scattering  in  the  detector  which 
had  an  energy  below  300  keV.  In  this  manner  the  detector  response  to  low 
energy  neutrons  was  estimated  in  the  calculations.  To  avoid  infinite  vari- 
ance results,  scattering  within  one  centimeter  of  the  detector  center  scored 

2 

with  the  average  value  of  1/r  over  the  sphere  (in  this  case  3)  rather  than 
2 

with  1/r  which  could  be  very  large.  The  random  walk  of  the  forward  17 
meter  calculation  was  maintained  so  that  the  effect  of  capture  in  the  detector 
would  not  be  clouded  by  statistics.  Figure  9 compares  the  two  forward  17 
meter  steel/styrofoam  calculations  with  and  without  capture  in  the  detector 
with  the  experiment.  Notice  that  within  statistics  the  differences  between 
the  experiment  and  the  calculation  have  vanished,  the  calculation  and  ex- 
periment agree.  Thus,  it  appears  that  for  the  steel  styrofoam  sample, 
capture  in  the  detector  accounts  for  30  to  40  percent  of  the  late  time  dose. 
For  budgeting  reasons  it  was  not  possible  to  repeat  the  other  calculations, 
including  the  Sprint  missile  calculation,  but  it  appears  this  effect  will  ex- 
plain why  the  calculated  late  time  results  are  usually  lower  than  the  experi- 
ments. 

Because  the  detector  is  thin  to  all  energy  neutrons  it  is  possible  to 
estimate  the  dosimeter  low  energy  neutron  response  from  the  multigroup 
homogenized  detector  cross  sections.  The  necessary  equation  is  derived 
as  follows:  Because  the  detector  is  thin, 


is  the  number  of  gamma  rays  generated  in  the  detector  volume  assuming 
0.  is  the  neutron  flux  in  group  i at  the  detector.  A and  i are  the  mean 
projection  and  the  mean  chord  length  of  the  detector  respectively  so  that 
SJeV,  the  detector  volume.  £*  is  the  neutron  macroscopic  total  cross 
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section  in  units  cm  and  G.  is  gamma  ray  generation  probability  given  a 
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TIME  AFTER  LINAC  BURST  (MICROSECONDS) 

Figure  9.  Effect  of  Neutron  Capture  in  the  Dosimeter.  Comparison 
of  17  Meter  Steel  Styrofoam  Forward  Calculations  with 
and  without  Neutron  Capture  in  the  Dosimeter. 


neutron  event  (G.  is  GAMGEN  in  MORSE).  This  source  of  gamma  rays 
can  be  converted  to  a gamma  ray  flux  by  multiplying  by  the  average  track 
length  (in  this  case  i /2  because  the  gamma  source  is  inside  the  detector) 
divided  by  the  volume.  This  gamma  flux  is  converted  to  a dose  by  multiply- 
ing by  D =Z^f.  .D  where  f. . is  the  fraction  of  gamma  rays  generated  in 
1 j J !J 

group  j from  neutrons  in  group  i and  IT  is  the  dosimeter  response  to 
a gamma  ray  flux  in  group  j.  The  detector  response  to  low  energy  neu- 
trons thus  becomes 


where  the  mean  chord  length  7 is  given  by  4V/A  = 2RL/(R+L)  = 1.  77  cm 
(A  here  is  the  total  surface  area  of  the  detector).  Table  16  presents 
£■  > G.,  D.,  and  R.  for  all  low  energy  neutron  groups.  Notice  that  the 
epithermal  neutron  response  is  considerably  lower  than  the  fast  neutron 
response  which  was  presented  in  Table  5,  but  that  the  response  in  the 
bottom  three  groups  certainly  is  not  negligible. 
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TABLE  16.  CALCULATION  OF  THE  DETECTOR  RESPONSE 
TO  LOW  ENERGY  NEUTRONS 
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5.  CALCULATION  OF  THE  SPRINT  MISSILE  EXPERIMENTS 


In  this  chapter  the  real  purpose  of  the  program,  that  is  to  verify 
our  ability  to  calculate  internal  doses  in  complicated  three  dimensional 
missiles,  is  addressed.  To  gain  confidence  in  our  geometry  model  a pre- 
liminary experiment  on  the  autopilot  was  performed  with  the  dosimeter 
137 

in  place  and  a Cs  source  located  at  eight  positions  around  the  autopilot. 
The  first  section  of  this  chapter  discusses  an  adjoint  calculation  of  this 
experiment  which  shows  that  with  one  exception  the  geometry  model  of 
the  autopilot  is  as  accurate  as  was  hoped.  The  last  two  sections  of  the 
■ chapter  discusses  our  calculations  of  the  Sprint  missile  in  the  linac  source 

with  the  detector  in  the  MGS  and  the  autopilot  respectively.  Eight  measure- 
ments of  the  missile  were  performed  which  were  composed  of  two  detector 
positions  (autopilot  and  MGS),  two  flight  path  lengths  (53.095  meters  and 
• , 17.  59  meters),  and  two  missile  orientations  (90°  and  30°  off  the  nose). 

Because  of  budgetary  consideration  only  five  of  these  were  calculated; 
these  were  the  four  high  energy  experiments  at  53.095  meters,  and  only 
the  90°  MGS  experiment  at  17.  59  meters. 

1 

5.  1 CALCULATION  OF  THE  Cs  AUTOPILOT  EXPERIMENT 
137 

The  Cs  experiment  is  depicted  in  Fig.  10  which  was  supplied  by 

IRT.  The  response  of  dosimeter  inside  the  autopilot  to  a 0.662  MeV  iso- 

137 

tropic  gamma  ray  source  (Cs  ) was  measured  at  each  of  the  eight  posi- 
tions shown  in  the  figure.  The  experimental  results  have  »ractional  stan- 
dard deviations  less  than  0.4%. 
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Figure  10.  Configuration  of  the  Calculation  of  the  Sprint  A/P 
Response  to  137  Cs  Gamma  Rays  (Not  to  Scale). 
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We  were  able  to  calculate  the  experiment  with  a single  adjoint 

calculation  which  employed  a next  flight  estimator  to  each  of  the  eight 

source  positions.  The  calculation  used  twelve  groups  tailored  especially 
137 

for  Cs  sources.  Table  8 presents  the  group  structure  and  the  dos- 
imeter response.  The  calculation  tracked  2,000  histories  in  4 minutes 
on  the  CDC  7600  computer.  This  slow  run  time  resulted  from  eight  next 
flight  estimations  in  a 46  region  geometry.  The  calculation  generated 
7,  400  scatterings  and  estimated  the  scattered  component  to  10  percent 
fractional  standard  deviation.  The  uncollided  component  has  no  variance 
because  both  the  source  and  dosimeter  were  assumed  to  be  points,  and  the 
source  energy  was  ail  in  the  first  energy  group. 

Table  17  compares  our  calculation  to  the  experiments.  Notice  that 
except  for  detectors  2 and  3,  the  calculation  is  systematically  15  to  35  per- 
cent lower  than  the  experiment.  This  is  probably  because  our  homogenized 
model  allows  no  streaming  between  components  which  certainly  will  happen 
with  x-rays  incident  on  real  missile  hardware. 

Another  effect  that  could  cause  the  calculations  to  be  low  is  that 
137 

the  IRT,  Cs  source  was  held  in  a lead  shield  which  could  cause  an 
additional  component  to  the  source  from  lead  backscattered  gamma  rays. 
The  IRT  source  strength  included  only  the  667  keV  gamma  rays.  However, 
these  backscattered  gamma  rays  would  have  a considerably  lower  energy 
and  probably  could  not  significantly  penetrate  the  Autopilot.  The  large 
difference  between  the  experiment  and  calculation  for  the  second  source 
position  (-y  axis)  indicates  a mass  error  in  one  of  the  regions,  probably 
region  9 (the  Structural  Filter).  This  discrepancy  however  was  not  re- 
solved, and  the  geometry  model  as  presented  in  Chapter  3 was  used  for 
the  rest  of  the  missile  calculations. 
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TABLE  17.  COMPARISON  OF  Cs  137  CALCULATION 
WITH  EXPERIMENT 


Source 

Position 

Rad  (Si) 

IRT 

Experiment 

% Difference 
Calc-Exp 
Exp 

Uncollided** 

Scattered 

Total* 

1 

6.  71  (-16) 

5.76  (-16) 

12.47  (-16) 

10.95  (-16) 

+ 14 

2 

4.39 

5.22 

9.  61 

6.235 

+ 54 

3 

4.25 

5.23 

9.48 

4.090 

+ 131 

4 

4.33 

4.75 

9.08 

7.776 

+ 17 

5 

5.86 

5.34 

11.20 

10.78 

+4 

6 

3.94 

4. 10 

8.01 

9.061 

-11 

7 

3.51 

3.46 

6.97 

7.328 

-5 

8 

4.  63 

4.42 

9.05 

7.972 

+ 13 

‘Statistics  are  10-15%  on  the  scattered  contribution. 
“Uncollided  with  no  A/P  present  would  be  1.888  (-15) 
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137 

The  conclusions  of  the  Cs  experiments  are  that  within  a 10  to 
30  percent  range  caused  by  x-ray  streaming  between  components,  the 
autopilot  model  is  accurate  except  for  one  direction.  This  discrepancy 
could  be  corrected  with  the  addition  of  mass  and/or  high  atomic  number 
elements  to  the  model  in  region  9.  The  effect  of  streaming  between 
components  should  not  affect  our  ability  to  calculate  neutron  and  gamma 
ray  doses  as  this  radiation  will  easily  penetrate  individual  components. 

5.  2 CALCULATIONS  OF  THE  SPRINT  MISSILE  EXPERIMENT 

Most  of  the  details  pertinent  to  the  Sprint  missile  calculations 
have  already  been  presented.  The  calculational  technique  was  essentially 
the  adjoint  time  dependent  method  discussed  in  Section  4.2.  The  geom- 
etry model  was  documented  in  Section  3.  3 and  the  source  and  dosimeter 
treatments  were  discussed  in  Section  3.  4.  The  only  remaining  detail  is 
to  define  the  source  beam  orientations  for  each  of  the  calculations  which 
are  presented  in  Table  18.  The  adjoint  next  flight  estimation  to  the  beam 
assumed  that  the  beam  began  at  a point  at  the  location  to  the  Be  photoneu- 
tron target.  Rather  than  rotating  the  missile  for  each  calculation,  it  was 
much  simpler  to  rotate  the  beam.  Therefore,  Table  18  gives  the  coordi- 
nate of  the  source  point  and  the  z intercept  of  the  beam  center  wilh  the 
missile  coordinates  for  each  of  the  five  missile  calculations  we  performed. 

Table  19  presents  some  of  the  Monte  Carlo  details  of  these  calcula- 
tions which  were  all  performed  on  a CDC  7600  computer.  Notice  that  each 
of  the  runs  were  about  10  minutes  for  20,000  adjoint  histories.  Also 
notice  that  the  90°  and  30°  calculations  of  the  same  detector  used  exactly 
the  same  random  walk,  but  because  of  estimation  to  different  beam  orien- 
tations, results  and  run  times  are  different. 
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TABLE  18.  COORDINATES  OF  THE  SOURCE  BEAM 
FOR  THE  SPRINT  CALCULATION 


1 

Source  Coordinates  (cm) 

Z Intercept 
(cm) 

Calculation 

X 

y 

z 

53 

AP 

90° 

-5309.5 

0.0 

136.8 

136.8 

53 

MGS 

90° 

-5309.5 

0.0 

100.91 

100.91 

53 

AP 

CO 

o 

o 

-2654.7 

0.0 

-4457.6 

128.  5 

53 

MGS 

30° 

-2666.7 

0.0 

-4490.4 

128.  5 

17 

MGS 

co 

o 

o 

-1759.0 

0.0 

100.91 

100.91 
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TABLE  19.  MONTE  CARLO  DETAILS  OF  SPRENT  MISSILE  CALCULATION 
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The  time  dependent  results  of  these  calculations  are  compared  with 
the  IRT  experimental  results  in  figures  11  through  15  where  the  vertical 
error  bars  represent  one  standard  deviation  of  the  calculation.  Figures  11 
and  12  show  the  53  meter  MGS  results  at  the  90°  and  30°  orientations.  The 
agreement  shown  here  is  remarkable,  strongly  demonstrating  our  ability 
to  calculate  the  high  energy  neutron  and  secondary  gamma  ray  silicon 
doses  in  missile  electronics.  This  also  demonstrates  that  the  eleven 
region  MGS  model  is  certainly  an  accurate  enough  geometrical  representa- 
tion for  missile  n,  y calculations. 

Figure  13  shows  the  results  of  the  17  meter  MGS  calculation  at 
90°.  The  results  here  are  dominated  by  neutron  thermalization  in  the 
missile  generating  secondary  gamma  rays  from  neutron  capture.  The 
results  shown  are  consistently  low  in  the  same  manner  that  the  steel-styro- 
foam 17  meter  calculations  were  low.  We  have  already  discussed  that  this 
effect  is  due  to  using  a dosimeter  response  that  neglected  neutron  capture 
in  the  dosimeter. 

Figures  14  and  15  show  the  53  meter  autopilot  results  at  the  90° 
and  30°  orientations  respectively.  Here  there  appears  to  be  a systematic 
difference.  Before  2.8  m seconds  the  agreement  is  very  good,  but  from 
2.8  m seconds  to  10  p seconds  the  calculations  are  nearly  a factor  of  two 
higher  than  the  experiment.  The  cause  of  this  discrepancy  has  not  been 
determined,  but  it  certainly  could  be  caused  by  an  error  in  the  cross  section 
mixtures.  A preliminary  examination  of  the  cross  section  mixtures  showed 
that  an  order  of  magnitude  error  in  the  copper  density  in  region  45  had 
always  existed  in  the  model  and  we  hoped  that  it  was  causing  this  difference. 
However,  a recalculation  of  the  autopilot  at  53  meters  and  90°  showed 
nearly  the  same  results.  Therefore  at  this  time,  the  source  of  this 
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DOSE  RATE  (RAD  (Sil/MICROSECOND  - SOURCE  MONITOR) 


TIME  AFTER  UNAC  BURST  (MICROSECONDS) 


Figure  11.  Comparison  of  the  53  Meter  Missile 
Guidance  Set  Adjoint  Calculation  at 
90°  With  the  Linac  Experiment 
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Figure  12.  Comparison  of  the  53  Meter  Missile 
Guidance  Set  Adjoint  Calculation  at 
30°  With  the  Linac  Experiment 


TIME  AFTER  LINAC  BURST  (MICROSECONDS) 


Figure  13.  Comparison  of  the  17.  59  Meter  Missile  Guidance  Set 
Adjoint  Calculation  at  90°  With  the  Linac  Experiment. 
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TIME  AFTER  LINAC  BURST  (MICROSECONDS) 


Figure  14.  Comparison  of  the  53  Meter  Autopilot 
Adjoint  Calculation  at  90°  with  the 
Linac  Experiment 


DOSE  RATE  (RAO  (Si)/MICRO$ECOND  - SOURCE  MONITOR) 


TIME  AFTER  LINAC  BURST  (MICROSECONDS) 


Figure  15.  Comparison  of  the  53  Meter  Autopilot 
Adjoint  Calculation  at  30°  with  the 
Linac  Experiment 
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difference  is  not  known.  However  it  is  clear  that  this  difference  would 
only  make  a small  change  in  any  calculation  of  the  total  neutron  induced 
missile  dose  which  is  dominated  by  both  high  and  low  energy  neutrons. 
Neutrons  in  the  range  from  10  keV  to  1 MeV  which  are  causing  this  dif- 
ference have  the  least  affect  on  total  dose. 

Several  conclusions  can  be  drawn  from  these  calculations.  First, 
the  silicon  dose  for  high  energy  neutrons  is  accurately  calculated  in  all 
cases.  Second,  the  eleven  region  MGS  model  gives  better  agreement  with 
experiment  than  the  46  region  autopilot  model.  This  indicates  that  the 
eleven  regions  are  certainly  adequate  and  that  perhaps  by  going  to  greater 
detail  we  have  introduced  errors  in  the  model.  Third,  both  the  MGS  and 
autopilot  results  are  similar,  both  in  magnitude  and  shape  of  the  time 
dependent  results.  This  could  mean  that  the  effects  of  the  various  elements 
average  out  and  that  in  fact  reasonable  estimates  of  internal  doses  could  be 
made  without  as  detailed  geometrical  modeling.  These  conclusions  will  be 
discussed  further  in  the  final  chapter. 
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6.  CONCLUSIONS 


Overall  this  joint  experimental-calculational  program  has  been 
successful  in  demonstrating  that  missile  doses  generated  by  incident  neu- 
trons including  the  effects  of  inelastic  and  capture  gamma  rays  can  be  ac- 
curately calculated  with  existing  Monte  Carlo  techniques,  cross  sections 
and  missile  models.  Also  many  of  the  trends  shown  in. the  program  in- 
dicate that  these  same  doses  may  be  reasonably  estimated  with  simpler, 
less  expensive  techniques.  For  example  the  similarity  of  the  results  from 
the  MGS,  the  autopilot  and  the  steel-styrofoam  sample  indicate  that  general 
sizes  and  average  materials  dominate  the  results  rather  than  individual 
immediate  components  which  are  very  different  for  each  calculation  or  mea- 
surement. 

The  initial  stages  of  this  program  involved  designing  the  steel-styro- 
foam sample,  processing  cross  sections  into  multigroup  structure  and  pre- 
paring the  66  region  of  the  Sprint  missile.  The  multigroup  cross  section 
processing  proved  particularly  difficult.  As  we  attempted  to  process  new 
cross  section  evaluations,  modifications  had  to  be  made  to  our  version  of 
AMPX  to  correct  errors  which  had  not  been  previously  tested.  Newer  cross 
section  evaluations  caused  continual  updating  of  our  multigroup  library 
especially  as  the  newer  DNA  evaluations  become  available.  This  constant 
effort  was  necessary  in  order  to  assure  that  no  differences  between  experi- 
ment and  calculation  would  result  from  an  outdated  cross  section  set  or  a 
failure  to  use  an  available  gamma-ray  generation  set. 
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The  generation  of  the  66  region  sprint  model  was  greatly  simplified 
by  the  outstanding  job  performed  by  Martin  Marietta  giving  us  their  present 
model  and  documenting  in  great  detail  the  regions  where  the  experimental 
and  flight  configurations  of  missile  differed,  particularly  where  electronics 
had  to  be  drilled  out  to  allow  insertion  of  the  dosimeter.  The  fact  that  this 
program  demonstrates  that  reasonable  estimates  are  possible  with  coarser 
modeling  could  save  funds  on  future  assessments  of  internal  missile  doses. 

The  steel-styrofoam  sample  experiment  has  proved  to  be  an  excel- 
lent benchmark  experiment.  Good  agreement  between  a time  dependent 
experiment  and  two  independent  calculational  techniques  at  all  times  serves 
to  verify  not  only  experimental  techniques  but  also  the  cross  sections  and 
the  calculation  methods.  Previous  comparisons  of  forward  and  adjoint 
Monte  Carlo  calculations  have  verified  that  the  two  methods  generate  the 
same  total  dose,  but  this  single  number  often  is  dominated  by  only  one 
aspect  of  the  transport  problem.  Agreement  from  early  to  late  times  of 
the  forward  and  adjoint  calculations  adds  considerable  confidence  to  the 
relatively  new  adjoint  Monte  Carlo  method. 

Our  technique  of  performing  time  dependent  adjoint  calculations 
discussed  in  Section  4. 2 should  be  noted,  as  the  inclusion  of  time  and 
energy  dependent  sources  in  an  adjoint  calculation  is  not  straightforward 
and  involves  some  rather  tricky  code  modifications.  In  forward  Monte 
Carlo  calculations  time  and  energy  dependent  sources  are  generally  handled 
by  random  sampling  both  the  energy  and  time  from  the  appropriate  dis- 
tribution. However,  in  adjoint  Monte  Carlo,  scoring  a time  and  energy 
dependent  source  requires  a careful  mapping  of  the  source  time  and  the 
scoring  times  with  the  source  energy.  Our  discussion  in  Section  4. 2 
should  aid  other  code  users  attempting  similar  problems. 
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The  late  time  response  of  the  steel-styrofoam  cylinder  (Figure  8) 
shows  that  the  calculations  are  30-50%  lower  than  the  experiment.  How- 
ever, inclusion  in  the  calculation  of  neutron  capture  in  the  detector 
(Figure  9)  shows  excellent  agreement  with  the  experiment.  Because  the 
remainder  of  the  calculations  likewise  do  not  include  a dosimeter  response 
to  low  energy  neutrons,  they  also  should  be  somewhat  low  at  late  time. 

This  effect  is  seen  in  the  MGS  calculation  at  17  meters  (Figure  13). 
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Calculation  of  the  eight  Cs  measurements  on  the  autopilot  show 
that  for  seven  out  of  the  eight  orientations,  the  autopilot  model  appears 
accurate.  The  calculations  are  somewhat  lower  than  the  experiment,  but 
this  can  be  attributed  to  two  effects.  First,  the  autopilot  model  is  a series 
of  homogenizations  having  no  streaming  paths  between  them.  In  reality 
the  autopilot  contains  hundreds  of  discrete  components  surrounded  by  low 
density  foam  allowing  x-rays  to  scatter  around  the  components  avoiding 
absorption.  In  this  manner  a homogenized  model  will  give  underestimates 
of  the  dose.  Second,  the  calculation  did  not  include  any  low  energy  gamma 
rays  that  may  have  been  backscattered  from  the  lead  source  holder.  This 
effect  is  probably  negligible.  In  one  direction  the  Cs  experiment  indi- 
cates a mass  difference  in  the  autopilot,  but  this  discrepancy  was  not 
resolved  and  may  be  relative  to  the  differences  between  the  autopilot  and 
Linac  comparisons. 

The  actual  missile  calculations  showed  excellent  agreement  at  the 
MGS  for  all  energy  neutrons  and  for  high  energy  neutrons  in  the  autopilot. 
For  neutrons  from  1. 8 to  . 15  MeV  in  the  autopilot,  calculations  were  higher 
than  the  experiment.  This  probably  was  due  to  inaccuracies  in  the 
geometry  model.  The  good  results  with  the  eleven  region  MGS,  and  the 
discrepancies  with  the  46  region  autopilot  model  seems  to  indicate  that 
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increased  modeling  effort  does  not  generate  more  accurate  dose  calcula- 
tions, and  in  fact,  increases  the  changes  of  an  error.  In  spite  of  the  dif- 
ference at  the  autopilot  it  must  be  emphasized  that  calculations  of  the 
total  dose  would  still  be  accurate  because  neutrons  in  the  range  of  the 
differences  are  the  least  important  for  calculating  total  dose. 

Two  final  conclusions  about  future  missile  internal  dose  calculations 
can  be  made.  First  the  problem  we  found  with  the  dosimeter  response 
indicates  that  care  should  be  taken  in  the  immediate  region  of  a com- 
ponent whose  dose  is  to  be  calculated.  The  best  approach  would  be  a 
preliminary  calculation  of  a component  response  that  includes  its  immediate 
region.  This  component  response  could  then  be  used  in  a coarse  missile 
model.  This  approach  would  include  the  detail  where  required,  yet  allow 
a single  missile  model  for  fast  and  accurate  calculations.  The  component 
response  might  include  a volume  of  material  whose  mean  chord  length 
would  be  approximately  . 2 mean  free  paths  of  2 Me V neutrons.  This  is 
small  enough  that  single  scattering  would  dominate,  yet  large  enough  that 
capture  in  the  immediate  area  would  be  accurately  determined.  The 
component  response  could  be  calculated  by  adjoint  Monte  Carlo  with  a 
model  of  the  immediate  region. 

Finally  the  similarity  of  all  the  measurements  and  calculations  of 
this  program  (steel-styrofoam,  MGS  and  autopilot)  indicate  that  the  internal 
doses  are  governed  mostly  by  average  effects  such  as  an  average  homo- 
genization, average  size,  and  probably  average  hydrogen  content.  Then  it 
should  be  possible  to  generate  reasonably  accurate  internal  neutron  and 
gamma  ray  doses  from  the  gross  parameters  of  the  missile  without  either 
generating  detailed  missile  model  or  performing  expensive  adjoint  Monte 
Carlo  calculations. 
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